The gravitational wave burst signal from core collapse of rotating stars 
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We present results from detailed general relativistic simulations of stellar core collapse to a proto- 
neutron star, using two different microphysical nonzero-temperature nuclear equations of state as 
well as an approximate description of deleptonization during the collapse phase. Investigating a 
wide variety of rotation rates and profiles as well as masses of the progenitor stars and both equa- 
tions of state, we confirm in this very general setup the recent finding that a generic gravitational 
wave burst signal is associated with core bounce, already known as type I in the literature. The 
previously suggested type II (or "multiple-bounce") waveform morphology does not occur. Despite 
this reduction to a single waveform type, we demonstrate that it is still possible to constrain the pro- 
genitor and postbounce rotation based on a combination of the maximum signal amplitude and the 
peak frequency of the emitted gravitational wave burst. Our models include to sufficient accuracy 
the currently known necessary physics for the collapse and bounce phase of core-collapse super- 
novae, yielding accurate and reliable gravitational wave signal templates for gravitational wave data 
analysis. In addition, we assess the possiblity of nonaxisymmetric instabilities in rotating nascent 
proto-neutron stars. We find strong evidence that in an iron core-collapse event the postbounce 
core cannot reach sufficiently rapid rotation to become subject to a classical bar-mode instability. 
However, many of our postbounce core models exhibit sufficiently rapid and differential rotation to 
become subject to the recently discovered dynamical instability at low rotation rates. 

PACS numbers: 04.25.D-, 04.30.Db, 97.60.Bw, 02.70.Bf, 02.70.Hm 



I. INTRODUCTION 

The final event in the life of a massive star is the catas- 
trophic collapse of its central, electron-degenerate core 
composed of iron-peak nuclei. When silicon shell burn- 
ing pushes the iron core over its effective Chandrasekhar 
mass, collapse is initiated by a combination of electron 
capture and photo-disintegration of heavy nuclei, both 
leading to a depletion of central pressure support. Mas- 
sive stars in the approximate mass range of about 10 
to 100 solar masses (M Q ) experience such a collapse 
phase until their homologously contracting [1, 2] inner 
core reaches densities near and above nuclear saturation 
density where the nuclear equation of state (EoS) stiff- 
ens, leading to an almost instantaneous rebound of the 
inner core (core bounce) into the still supersonically in- 
falling outer core. The hydrodynamic supernova shock is 
born, travels outward in radius and mass, but rapidly 
loses its kinetic energy to the dissociation of infalling 
iron-group nuclei and to neutrinos that deleptonize the 
immediate postshock material and stream off from these 
regions quasi-frecly. The shock stalls, turns into an ac- 
cretion shock and must be revived to produce the observ- 
able explosion associated with a core-collapse supernova. 
Mechanisms of shock revival are still under debate (a re- 
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cent review is presented in [3], but see also [4-6]) and 
may involve heating of the postshock region by neutri- 
nos, multi-dimensional hydrodynamic instabilities of the 
accretion shock, in the postshock region, and/or in the 
proto-neutron star, rotation, magnetic fields, and nuclear 
burning. If the shock is not revived, black-hole formation 
(on a timescale of ~ 1-2 s [7] ) is inevitable and the stel- 
lar collapse event may remain undetected by conventional 
astronomy or, perhaps, appear as a gamma-ray burst if 
the progenitor star has a compact enough envelope and 
sufficiently rapid rotation in its central regions [8, 9]. 

Conventional astronomy can constrain core-collapse 
supernova theory and the supernova explosion mecha- 
nism via secondary observables only, e.g., the explo- 
sion energy, ejecta morphology, nucleosynthesis yields, 
residue neutron star or black hole mass and proper mo- 
tion, and pulsar magnetic fields. Neutrinos and gravita- 
tional waves, on the other hand, are emitted deep inside 
the supernova core and travel to observers on Earth prac- 
tically unscathed by intervening material. They can act 
as messengers to provide first-hand and live dynamical 
information on the intricate multi-dimensional dynam- 
ics of the proto-neutron star and postshock region and 
may constrain directly the core-collapse supernova mech- 
anism. Importantly, core-collapse events that do not pro- 
duce the canonical observational astronomical signature 
or whose observational display is shrouded from view can 
still be observed in neutrinos and gravitational waves if 
occurring sufficiently close to Earth. 

Gravitational waves, in contrast to neutrinos, have not 
yet been observed directly, but an international array of 
gravitational wave observatories (see, e.g., [10]) is active 
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and taking data. Since gravitational waves from astro- 
physical sources are expected to be weak, their detec- 
tion is notoriously difficult and involves extensive signal 
processing and detailed analysis of the detector output. 
Chances for the detection of an astrophysical event of 
gravitational wave emission are significantly enhanced if 
accurate theoretical knowledge of the expected gravita- 
tional wave signature from such an event is at hand. 

Theoretical predictions of the gravitational wave sig- 
nature from a core-collapse supernova are complicated, 
since the emission mechanisms are very diverse. While 
the prospective gravitational wave burst signal from the 
collapse, bounce, and the very early postbounce phase 
is present only when the core rotates [11-18], gravita- 
tional wave signals with sizeable amplitudes can also be 
expected from convective motions at postbounce times, 
instabilities of the standing accretion shock, anisotropic 
neutrino emission, excitation of various oscillations in the 
proto-neutron star, or nonaxisymmetric rotational insta- 
bilities [17, 19-23]. 

In the observational search for gravitational waves 
from merging black hole or neutron star binaries, power- 
ful data analysis algorithms such as matched filtering are 
applicable, as the waveform from the inspiral phase can 
be modeled with high accuracy (see, e.g., [24]) and gravi- 
tational wave data analysts already have access to robust 
template waveforms that depend only on a limited num- 
ber of macroscopic parameters. In contrast, the complete 
gravitational wave signature of a core-collapse supernova 
cannot be predicted with template-level accuracy as the 
postbounce dynamics involve chaotic processes (turbu- 
lence, [magneto-] hydrodynamic instabilities) that are 
sensitive not only to a multitude of precollapse parame- 
ters, but also to small-scale perturbations of any of the 
hydrodynamic variables. 

While the complete supernova gravitational wave sig- 
nature may remain inaccessible to template-based data 
analysis, a number of individual constituent emission pro- 
cesses, in particular those involving coherent global bulk 
dynamics and/or rotation, allow, in principle, for accu- 
rate and robust waveform predictions that may be ap- 
plied to template-based searches. Rotating core collapse 
and core bounce as well as pulsations or nonaxisymmet- 
ric rotational deformations of a proto-neutron star con- 
stitute this group of processes. Among them, rotating 
collapse and bounce is the historically most extensively 
studied case (see, e.g., [25] for a historical review) and 
may be the most promising for becoming robustly pre- 
dictable in its gravitational wave emission. Yet, to date, 
the gravitational wave signal from rotating stellar core 
collapse and bounce has not been predicted with the de- 
sired accuracy and robustness. 

These deficiencies of previous simulations result from 
the fact that the physically realistic modeling of core col- 
lapse requires a general relativistic description of consis- 
tently coupled gravity and hydrodynamics in conjunction 
with a microphysical treatment of the sub- and supernu- 
clcar EoS, electron capture on heavy nuclei and free pro- 



tons, and neutrino radiation transport. Only very few 
multi-dimensional general relativistic codes have recently 
begun to approach these requirements [17, 18]. In addi- 
tion, the properties of the EoS around and above nuclear 
density are not very well constrained by theory or experi- 
ments. The same applies to the rotation rate and angular 
velocity profile of the progenitor core, which are also not 
directly accessible by observation and very difficult to 
model numerically in stellar evolution codes. Further- 
more, variations with progenitor structure and mass are 
to be expected. Therefore, the influence of rotation and 
progenitor structure on the collapse and bounce dynam- 
ics and thus the gravitational wave burst signal must be 
investigated by extensive and computationally expensive 
parameter studies. 

Previous parameter studies have considered a large va- 
riety of rotation rates and progenitor core configurations, 
but generally ignored important microphysical aspects 
and/or the influence of general relativity. Monchmeyer 
ct al. [12] performed axisymmctric Newtonian calcula- 
tions with progenitor models from stellar evolutionary 
studies. They employed the microphysical nuclear EoS 
of Hillcbrandt and Wolff [26] and included deleptoniza- 
tion via a neutrino leakage scheme and electron cap- 
ture on free protons. Capture on heavy nuclei was ne- 
glected, which resulted in a too high electron fraction 
Y e at core bounce and a consequently overestimated in- 
ner core mass [2, 27]. In that study a limited set of 
four calculations was computed and two qualitatively and 
quantitatively different types of gravitational wave burst 
signals were identified. Their morphology can be clas- 
sified alongside with the collapse and bounce dynamics: 
Type I signals are emitted when the collapse of the quasi- 
homologously contracting inner core is not strongly influ- 
enced by rotation, but stopped by a pressure-dominated 
bounce due to the stiffening of the EoS near nuclear den- 
sity p nuc w 2 x 10 14 g cm -3 , where the adiabatic in- 
dex 7 GOS rises above 4/3. This leads to a bounce with 
a maximum core density p max > p nU c- Type II signals 
occur when centrifugal forces, which grow during con- 
traction owing to angular momentum conservation, are 
sufficiently strong to halt the collapse, resulting in con- 
secutive (typically multiple) centrifugal bounces with in- 
termediate coherent re-expansion of the inner core, seen 
as density drops by sometimes more than an order of 
magnitude; thus here p ma x < Pnuc after bounce. Type I 
and II dynamics and waveforms were also found in the 
more recent Newtonian studies by Kotake et al. [15], who 
employed a more complete leakage/capture scheme, but 
still obtained too high Y e at bounce, and by Ott et al. [16], 
who performed an extensive parameter study and for the 
first time also considered variations in progenitor star 
structure, but neglected dclcptonization during collapse. 

Zwerger and Miiller [13] carried out an extensive two- 
dimensional Newtonian study of rotating collapse of ide- 
alized polytropcs in rotational equilibrium [28] with a 
simplified hybrid EoS, consisting of a polytropic and a 
thermal component [29]. Electron capture during col- 
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lapse was mimicked by an instantaneous lowering of the 
adiabatic index 7 oos from its initial value of 4/3 to trig- 
ger the onset of collapse. At p nuc , the adiabatic index 
was raised to > 2 to qualitatively model the stiffening of 
the nuclear EoS. Zwerger and Miiller also obtained the 
previously suggested signal types and introduced type III 
signals that appear in a pressure-dominated bounce when 
the inner core has a very small mass due to very efficient 
electron capture (approximated in [13] via a 7 GOS < 1.29 
in their hybrid EoS). Obergaulinger et al. [30] also em- 
ployed the hybrid EoS, but included magnetic fields. 
They introduced the additional dynamics/signal type IV, 
which occurs only in the case of very strong precollapse 
core magnetization. They found that weak to moderate 
core magnetization in agreement with predictions from 
stellar evolution theory (see, e.g., [31]) has little effect 
on the collapse and bounce dynamics and the resulting 
gravitational wave signal. This finding is in agreement 
with [32] (see also [5, 33]), where magneto-rotational col- 
lapse simulations were performed, a smaller model set 
was considered, but the neutrino leakage scheme of [15] 
was employed and it was made use of two different mi- 
crophysical EoSs to study the EoS dependence of the 
collapse dynamics and gravitational wave signal. 

The first extensive set of general relativistic simula- 
tions of rotating iron core collapse to a proto-neutron 
star were presented by Dimmclmcicr ct al. [14], who em- 
ployed an analytic hybrid EoS and polytropic precollapse 
models in rotational equilibrium as initial data (but see 
also the pioneering early work of [34]). These simula- 
tions were subsequently confirmed in [25, 35-37]. Dim- 
melmeier et al. studied a subset of the models in [13] 
in the same parameter space of rotation rate and degree 
of differential rotation, and found that general relativistic 
effects counteract centrifugal support and shift the occur- 
rence of type II dynamics and wave signals to a higher 
precollapse rotation rate at a fixed degree of differential 
rotation. 

Recently, new general relativistic simulations of rotat- 
ing core collapse in two and three dimensions were carried 
out by Ott et al. [17, 25, 38] who included the microphys- 
ical EoS of Shen et al. [39] , precollapse models from stel- 
lar evolutionary calculations as well as an approximate 
deleptonization scheme [40]. The results of these calcu- 
lations indicate that the gravitational wave burst signal 
associated with rotating core collapse is exclusively of 
type I. In addition, the simulations showed that rotating 
stellar iron cores stay axisymmetric throughout collapse 
and bounce, and only at postbounce times develop non- 
axisymmetric features. 

In a general relativistic two-dimensional follow-up 
study, Dimmclmeier et al. [18, 41] considerably extended 
the number of models and comprehensively explored a 
wide parameter space of precollapse rotational configu- 
rations. Even for this more general setup they found 
gravitational wave signals solely of type I form, although 
for rapid precollapse rotation some of their models ex- 
perience a core bounce due to centrifugal forces only, 



which however is always a single centrifugal bounce rather 
than the multiple ones observed in earlier work (see, e.g., 
[13, 14, 16]). They identified the physical conditions that 
lead to the emergence of this generic gravitational wave 
signal type and quantified their relative influence. These 
results strongly suggest that the waveform of the gravi- 
tational wave burst signal from the collapse of rotating 
iron cores in a core-collapse event is much more generic 
than previously anticipated. 

In this work, we extend the above study of the gravi- 
tational wave signal from rotating core collapse and con- 
sider not only variations in the precollapse rotational con- 
figuration, but also in progenitor structure and nuclear 
EoS. In this way, we carry out the to-date largest and 
most complete parameter study of rotating stellar core 
collapse that includes all the (known) necessary physics 
to produce reliable predictions of the gravitational wave 
signal associated with rotating collapse and bounce. All 
our computed gravitational wave signals are made avail- 
able to the detector data analysis community in a freely 
accessible waveform catalog [42] . 

We perform a large number of two-dimensional sim- 
ulations with our general relativistic core-collapse code 
CoCoNuT and employ 11.2, 15.0, 20.0, and 40.0 M Q 
(masses at zero-age main sequence) precollapse stellar 
models from the stellar evolutionary studies of Heger et 
al. [31, 43]. In addition to the EoS by Shen et al. [39] 
used in our previous studies, we also calculate models 
with the EoS by Lattimer and Swesty [44]. We describe 
in detail and explain comprehensively the qualitative and 
quantitative aspects of the collapse and bounce dynamics 
and the resultant gravitational wave signal. We lay out 
the individual effects of general relativity, deleptoniza- 
tion, precollapse stellar structure and rotational configu- 
ration, and nuclear EoS on the gravitational wave signa- 
ture from rotating core collapse. We study the prospects 
for nonaxisymmetric rotational instabilities in our post- 
bounce cores, which could lead to an enhancement of the 
gravitational wave signature. Furthermore, we set our 
model gravitational radiation waveforms in context with 
present and future detector technology and assess their 
detectability. 

This paper is organized as follows: In Section II we 
introduce our treatment of the general relativistic space- 
time curvature and hydrodynamics equations. Further- 
more, we introduce our variants of the two microphysi- 
cal EoS we employ, the scheme for deleptonization and 
neutrino pressure contributions, our precollapse model 
set, and the gravitational wave extraction technique cm- 
ployed. Section III discusses the numerical methods used 
in the CoCoNuT code and the computational grid setup 
for the simulations presented in this paper. In Section IV 
we present the collapse dynamics and waveform morphol- 
ogy of our simulated models, while in Section V we in- 
vestigate the stratification of the postbounce core and 
its impact on the gravitational wave signal. The de- 
tection prospects for the gravitational wave burst from 
core bounce are discussed in Section VI, while the ro- 
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tational configuration of the proto-neutron star and its 
susceptibility to nonaxisymmetric rotational instabilities 
are examined in Section VII. Finally, in Section VIII, we 
summarize and discuss our results. 

Throughout the paper we use a spacelike signature 
(—,+,+,+) and units in which c = G = 1. Greek in- 
dices run from to 3, Latin indices from 1 to 3, and we 
adopt the standard Einstein summation convention. 



II. PHYSICAL MODEL AND EQUATIONS 
A. General relativistic hydrodynamics 



Here v 
9il V and 7y 



v l — P l /a, and g and 7 are the determinant of 
respectively, with y/—g = a^/7- T* v are the 
four-Christoffcl symbols. Since we use a microphysical 
EoS that requires information on the local electron frac- 
tion per baryon Y e , we add an advection equation for the 
quantity DY e to the standard form of the conservation 
equations (4). The radiation stress due to the neutrino 
pressure P v (as defined in Section II D), is included in 
the form of an additive term in the source of both the 
momentum and energy equations. Note also that here we 
use an analytically equivalent reformulation of the energy 
source term in contrast to the one presented in [14]. 



We adopt the Arnowitt-Deser-Misner (ADM) 3 + 1 
formalism of general relativity to foliate the spacetime 
endowed with a four-metric g^ into spacelike hypersur- 
faces [45]. In this approach the line element reads 



ds l 



-a 2 dt 2 + 7 y {dx l + (i l dt)(dx 3 + /3 j dt), (1) 



where a is the lapse function, (3 l is the shift vector, and 
7,j is the spatial three-metric induced in each hypersur- 
face. 

The hydrodynamic evolution of a perfect fluid in gen- 
eral relativity with four- velocity u^, rest-mass current 
J M = pu^, and stress-energy tensor T^ v = phu fJ, u lJ + 
Pg^ v is determined by a system of local conservation 
equations, 



0, 



0, 



(2) 



where V M denotes the covariant derivative with respect 
to the four-metric. Here p is the rest-mass density, 
h = 1 + e+P/p is the specific enthalpy, P is the fluid pres- 
sure, and the three-velocity with respect to an Eulerian 
observer moving orthogonally to the spacelike hypersur- 
faces is given by v % = u l /(au°) + (3 l /a. We define a set 
of conserved variables as 

D = pW, S l = phW 2 v\ t = phW 2 - P - D. (3) 



au° is the Lorentz factor, 

1/VT 



ViV 



In the above expressions W = 
which satisfies the relation W 

The local conservation laws (2) are written as a 
first-order, flux-conservative system of hyperbolic equa- 
tions [46], 



dt 



dx 1 



with 



U=[D,S jt T,DY e ], 



[Dv^Sjtf + 8)P, Tv l + Pv\DY e v l ] 



(4) 
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S = 



(J. 



d 9iiv dP v 



da 



2 ft* -F'^-^N 



B. Metric equations in the conformal flatness 
approximation 

Using the ADM 3+1 formalism, the Einstein equations 
split into a coupled set of first-order evolution equations 
for the three- metric 7^ and the extrinsic curvature Kij, 



d t j i:j = 2n A ,, + V t f3j + ViA, 
d t K i:j = -ViV.-Q + a (Rij - 2/\ ,/, A'*") + 
f3 k V k K l3 + K lk V 3 p k + K jk ViP k 
8-n-a (^S i:j - -y- (£* - padm) 

and constraint equations, 

= R-K ij K i 3 -16TTPADM, 



(8) 
(9) 



= ViAT' 



8ttS j . 



(10) 
(11) 



In the above equations, Vj is the covariant derivative with 
respect to the three-metric 7^ , Rij is the three- Ricci ten- 
sor and R is the scalar three-curvature. The projection of 
the stress-energy tensor onto the spatial hypersurface is 
Sij = phW 2 ViVj+jijP, the ADM energy density is given 
by Padm = phW 2 — P, and = phW 2 v L is the momen- 
tum density. In addition, we have chosen the maximal 
slicing condition for which the trace of the extrinsic cur- 
vature vanishes: K = 0. 

In order to simplify the ADM metric equations and 
to ameliorate the stability properties when numerically 
solving those equations, we employ the conformal flat- 
ness condition (CFC) introduced in [47] and first used in 
a pseudo-evolutionary context in [48]. In this approxi- 
mation the spatial three-metric is replaced by the con- 
formally flat three-metric, 7^ = <f> 7y, where 7^ is the 
flat-space metric and <p is the conformal factor. Then the 
metric equations (8-11) reduce to a set of elliptic equa- 
tions for 0, a, and /?*, 



L{a<j>) = 2na<b 5 ( E + 2S + 

l07T 



(12) 
(13) 



A/3 1 = 167ra0 4 S l + 2<j> 1Q K ij V~ - -V l V fe /3 fc , (14) 

0° 3 
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where A and V are the Laplace and covariant derivative 
operators associated with the flat three-metric, and S = 
j^Sij. The CFC metric equations (12-14) do not contain 
explicit time derivatives, and thus the metric components 
are evaluated in a fully constrained approach. 

Imposing CFC in a spherically symmetric spacetime is 
equivalent to solving the exact Einstein equations. For 
nonsphcrical configurations the CFC approximation may 
be roughly regarded as full general relativity without the 
dynamical degrees of freedom of the gravitational field 
that correspond to the gravitational wave content [49]. 
However, even spacetimes that do not contain gravita- 
tional waves can be not conformally flat. A prime exam- 
ple are the spacetime of a Kerr black hole [50] or rotat- 
ing fluids in equilibrium. For rapidly rotating models of 
stationary neutron stars the deviation of certain metric 
components from conformal flatness has been shown to 
reach up to ~ 5% in extreme cases [51], while the oscilla- 
tion frequencies of such models typically deviate even less 
from the corresponding values obtained in full general rel- 
ativistic simulations [52] . In the context of rotating stel- 
lar core collapse the excellent quality of the CFC approx- 
imation has been demonstrated extensively [17, 35, 36]. 

Due to its fully constrained nature, the CFC approxi- 
mation permits a straightforward and numerically more 
robust implementation of the metric equations in coor- 
dinate systems containing coordinate singularities (e.g., 
spherical polar coordinates) compared to a Cauchy free- 
evolution scheme. Furthermore, by definition it allows 
no constraint violations, which is a significant benefit in 
cases where a perturbation is added to the initial data. 
More details on the CFC equations can be found in, e.g., 
[14]. 



C. Equations of state 

In our simulations we employ two tabulated nonzero- 
temperature equations of state, the one by Shen et 
al. [39, 53] (Shen EoS), and the one by Lattimer and 
Swesty [44] (LS EoS). The LS EoS is based on a com- 
pressible liquid-drop model [54] . The transition from in- 
homogeneous to homogeneous matter is established by a 
Maxwell construction, and the nucleon-nucleon interac- 
tions are expressed by a Skyrme force. In our version 
of this EoS, the incompressibility modulus of bulk nu- 
clear matter is taken to be 180 McV and the symmetry 
energy parameter has a value of 29.3 MeV. In contrast, 
the Shen EoS is based on a rclativistic mean field model 
and is extended with the Thomas-Fermi approximation 
to describe the homogeneous phase of matter as well as 
the inhomogeneous matter composition. The parameter 
for the incompressibility of nuclear matter is 281 MeV 
and the symmetry energy has a value of 36.9 MeV. 

Both EoSs employed in this study are the same as in 
Marck ct al. [55] and include contributions of baryons, 
electrons, positrons, and photons. Furthermore, in this 
study the LS EoS has been extended to densities below 



p = 5.8 x 10 7 g cm -3 by a smooth transition to the Shen 
EoS which is tabulated down to p = 6.4 x 10 5 g cm~ 3 . 

The microphysical EoS returns the fluid pressure (and 
additional thermodynamic quantities) as a function of 
(p,T,Y e ), where T is the temperature. Since the hydro- 
dynamic equations (4) operate on the specific internal 
energy e, we determine the corresponding temperature T 
iteratively with a Newton-Raphson scheme and the EoS 
table. All interpolations are carried out in tri-linear fash- 
ion and the tables are sufficiently densely spaced to lead 
to an artificial entropy increase in an adiabatic collapse 
by not more than ~ 2%. 

D. Deleptonization and neutrino pressure 

Electron capture on free protons and heavy nuclei dur- 
ing collapse reduces Y e (i.e., "deleptonizes" the collapsing 
core) and consequently decreases the size of the homolo- 
gously collapsing inner core that depends on the average 
value of Y e in a roughly quadratic way (see, e.g., [56]). 
The material of the inner core is in sonic contact and 
determines the dynamics and the gravitational wave sig- 
nal at core bounce and in the early postbounce phases. 
Hence, deleptonization has a direct influence on the col- 
lapse dynamics and the gravitational wave signal, and 
thus it is essential to include deleptonization during col- 
lapse. 

Since multi-dimensional radiation-hydrodynamics cal- 
culations in general relativity are not yet computationally 
feasible, in our simulations we make use of a recently pro- 
posed approximative scheme [40] where deleptonization 
is parametrized based on data from detailed spherically 
symmetric calculations with Boltzmann neutrino trans- 
port, for which (as in [18]) we take the latest available 
electron capture rates [57]. Following the main assump- 
tion in [40] that the local electron fraction for each fluid 
element during the contraction phase can be modeled 
rather accurately by a dependence on the density only, 
these simulations yield a universal relation Y e (p). Fur- 
thermore, we find that this relation varies only slightly 
with progenitor mass, as shown in Fig. 1, where mod- 
els with identical progenitor but different EoS have the 
same color, but different hues (e.g., dark green versus 
light green for the s20 progenitor). Consequently, we uti- 
lize the 20.0 M Q progenitor to create such a profile Y e (p) 
for each of the two EoSs. This profile is then used to 
correct the value of Y e obtained from the advection by 
an amount 

AY e = min[0,F e (p) -Y e ] (15) 

after each time integration step. This procedure assures 
that Y e approaches the phenomenological input profile 
Y e (p) with the constraint that AY e must be negative. 
Accordingly, in order to model the entropy loss by neu- 
trinos escaping the collapsing core, for densities below an 
adopted neutrino trapping density p tT = 2x 10 12 g cm" 3 
the internal specific energy e is re-adjusted at constant 
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FIG. 1: Electron fraction Y e obtained from detailed spheri- 
cally symmetric calculations with Boltzmann neutrino trans- 
port versus the maximum density p max in the collapsing core. 
The EoS is encoded in dark hues for the Shen EoS and light 
hues for the LS EoS with the basis color specifying the pro- 
genitor mass. 



TABLE I: Properties of the iron core models used as initial 
data. Afp r0 g is the total zero-age main sequence mass of the 
progenitor star, M core and R CO ra are the mass and radius of 
the iron core, Mi and Ri are the mass and radius of the ini- 
tial model on the computational grid, and p c ,i is the prec- 
ollapse density at the center. The size of the iron core is 
determined by the condition that V e exceeds 0.497, while the 
initial model extends beyond the iron core to where the den- 
sity drops to 10 -4 p c ,i- pc,i deviates slightly from the original 
value of the models in [58] because of regridding to the more 
densely spaced central grid of the evolution code. 
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40.0 
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15.0 
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1.55 


1.83 
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1.56 


3.17 
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2.26 
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1.38 


1.35 


1.60 


3.18 


7.22 



p and Y e such that the specific entropy per baryon s is 
changed by 

As = -Ay / p ~ Mra + Me ~^ , (16) 

where E v = 10 McV is an average escape energy for the 
neutrinos, fce is the Boltzmann constant and where 
/i n , and /i e are the proton, neutron, and electron chemi- 
cal potentials, respectively. Note that when equilibrium 
between neutrinos and matter (i.e., /^-equilibrium) is es- 
tablished, this balance requires p\, v = p, p — fj, n + \i e for the 
neutrino chemical potential p, v . 

We stop deleptonization at the time of core bounce 
(i.e., as soon as the specific entropy s per baryon ex- 
ceeds 3&B at the outer boundary of the inner core). After 
core bounce, for lack of a simple yet accurate approxima- 
tion scheme for treating the further deleptonization in the 
nascent proto-neutron star, we advect Y e only passively 
according to the conservation equation (4), although this 
effectively prevents the factual cooling and contraction of 
the proto-neutron star. 

In all collapse phases, however, as in [40] we approxi- 
mate the pressure contribution of the neutrinos by that 
of an ideal Fermi gas, 

with i<3 being the Fermi-Dirac function of order 3. The 
neutrino pressure is included only in the regime which is 
optically thick to neutrinos, which we define for densities 
above pt T . 



E. Initial models 

All prcsupcrnova stellar models available to-date are 
end products of Newtonian spherically-symmetric stel- 
lar evolutionary calculations from hydrogen burning on 
the main sequence to the onset of core collapse by 
photo-dissociation of heavy nuclei and electron captures 
(see, e.g., [58]). Here, we employ various nonrotat- 
ing models of [58] with zero-age main sequence masses 
Mp r0 g = 11.2 Mq (core model sll. 2, here for simplic- 
ity labeled sll), 15.0 Af© (core model sl5), 20.0 Mq 
(core model s20), and 40.0 Mq (core model s40). Re- 
cently, the first presupernova models that include rota- 
tion in a one-dimensional approximate fashion have be- 
come available [31, 43], and of these we select ones with 
-^prog = 15.0 Mq (models el5a and el5b) as well as 
20.0 Mq (core models e20a and e20b). All progenitors 
have solar-mctallicity (at zero-age main sequence) and 
we generate our initial models by taking the data ob- 
tained from stellar evolution out to a radius i?i where 
the density drops to a value that equals 10 -4 of the ini- 
tial precollapse central density p c ^. Selected quantities 
that describe the properties of these stellar cores are sum- 
marized in Table I. 

We set those cores that are initially nonrotating (core 
models sll, sl5, s20, and s40) artificially into rotation 
according to the rotation law specified in [28] , where the 
specific angular momentum j is given by 

j=A 2 (Q c : 1 -Q). (18) 

Here the length A parametrizes the degree of differential 
rotation (stronger diffcrcntiality with decreasing A) and 
f2 Ci i is the precollapse value of the angular velocity at 
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the center. In the Newtonian limit, this reduces to 



Gravitational wave extraction 



— i 



A 2 



A 2 + t 2 sin 2 9 ' 



(19) 



with rsin^ being the distance to the rotation axis. 

In order to determine the influence of different angu- 
lar momentum distributions on the collapse dynamics, 
we parameterize the precollapse rotation of our models 
in terms of A (Al: A = 50,000 km, almost uniform; A2: 
A = 1, 000 km, moderately differential; A3: A = 500 km, 
strongly differential) and Q Ctl . The model nomenclature 
for the precollapse rotation parameters is shown in Ta- 
ble II. We have selected the rotational configuration of 
the models in such a way that for the s20 progenitor 
they are a representative subset of the models investi- 
gated in [18, 41]. They reflect different properties of the 
collapse dynamics and the gravitational radiation wave- 
form discussed in that work, namely pressure-dominated 
bounce with or without significant postbounce convective 
overturn as well as single centrifugal bounce. 

Note that models with the same rotation specification 
(but different progenitor mass or EoS) have an identical 
angular velocity profile, while the precollapse rotation 
rate /3; = Ti/|W|i, which is the precollapse ratio of ro- 
tational energy to gravitational energy, varies. We have 
decided to compare models with identical initial angu- 
lar velocity fi C) ; and not precollapse rotation rate as 
the latter quantity is rather sensitive to the chosen core 
radius R corc in the case of (almost) uniform rotation. 

The models that are based on progenitor calculations 
including rotation (core models el5a, el5b, e20a, and 
c20b) are mapped onto our computational grids under 
the assumption of constant rotation on cylindrical shells 
of constant distance to the rotation axis. We also point 
out that due to the one-dimensional nature, none of the 
considered models are in rotational equilibrium. Still, 
in slowly rotating initial models this effect is small and 
thus negligible. For more rapidly rotating models, which 
exhibits the strongest deviation from rotational equilib- 
rium, the collapse proceeds more slowly due to stabilizing 
centrifugal forces, and hence the star has more time for 
the adjustment to the appropriate angular stratifications 
for its rate of rotation. 

In this study, we focus on the collapse of massive 
presupernova iron cores with at most moderate differ- 
ential rotation and precollapse rotation rates that ex- 
cept for the most slowly rotating models lead to proto- 
neutron stars that are probably spinning too fast to yield 
spin periods of cold neutron stars in agreement with ob- 
servationally inferred injection periods of young pulsars 
into the P/P diagram [31, 59]. However, they may be 
highly relevant in the collapsar-type gamma-ray burst 
scenario [9, 59, 60]. 



We employ the Newtonian quadrupole formula in the 
first-moment of momentum density formulation as dis- 
cussed, e.g., in [14, 61, 62] to extract the gravitational 
waves generated by nonspherical accelerated fluid mo- 
tions. It yields the quadrupole wave amplitude A^ 2 as 
the lowest order term in a multipole expansion of the ra- 
diation field into pure-spin tensor harmonics [63]. The 
wave amplitude is related to the dimensionless gravita- 
tional wave strain h in the equatorial plane by 



h 



1^ ^20 



8.8524 x 10 



AE2 
-21 A 2Q 



10 kpc 



10 3 cm 



(20) 



where r is the distance to the emitting source. 

We point out that although the quadrupole formula is 
not gauge invariant and is only valid in the Newtonian 
slow- motion limit, for gravitational waves emitted by pul- 
sations of rotating NSs (i.e., in astrophysical situations 
comparable to collapsing stellar cores at bounce in terms 
of compactness and rotation rates) it yields results that 
agree very well in phase and to ~ 10-20% in amplitude 
with more sophisticated methods [61, 64]. 

In order to assess the prospects for detection by current 
and planned interferometer detectors and to specifically 
address the issue of detection range and expected event 
rates, we calculate the dimensionless characteristic grav- 
itational wave strain h c of the signal according to [65]. 
We first perform a Fourier transform of the gravitational 
wave strain h, 



h = 



s 2 * ift hdt. 



(21) 



To obtain the (detector dependent) integrated character- 
istic signal frequency 



fc 




(h 2 



Si, 



-df 



and the integrated characteristic signal strain 



'^{h 2 )fdf 

Oh 



1/2 



(22) 



(23) 



the power spectral density Sh of the detector is needed 
(with She = Sh(fc))- We approximate the average {h 2 } 
over randomly distributed angles by h , assuming op- 
timal orientation of the interferometer detector. From 
Eqs. (22, 23) the signal-to-noise ratio can be computed 
as SNR = /i c /[/&rms(/c)]j where h rms = fSh is the value 
of the rms strain noise (i.e., the theoretical sensitivity 
window) for the detector. 



III. NUMERICAL METHODS 

We perform all simulations using the CoCoNuT code 
described in detail in [14, 62]. The equations of general 
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TABLE II: Precollapse rotation properties of the core collapse models. A is the differential rotation length scale, S7 Cj i is the 
precollapse angular velocity at the center, and $ is the precollapse rotation rate. Note that the models el5a, el5b, e20a, 
and e20b have a rotation profile from the corresponding stellar evolution calculations, while onto all other models an artificial 
rotation profile is imposed. 
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relativistic hydrodynamics are solved in semi-discrete 
fashion. The spatial discretization is performed by means 
of a high-resolution shock-capturing (HRSC) scheme em- 
ploying a second-order accurate finite- volume discretiza- 
tion. We make use of the Harten-Lax-van Lecr-Einfeldt 
(HLLE) flux formula for the local Riemann problems and 
piecewise-parabolic reconstruction (PPM) of the primi- 
tive variables (p, u l ,e) at cell interfaces. For a review 
of such methods in general relativistic hydrodynamics, 
see [66]. The time integration and coupling with curva- 
ture are carried out with the method of lines [67] in com- 
bination with a second-order accurate explicit Rungc- 
Kutta scheme. Once the state vector U is updated in 
time, the primitive variables are recovered from the con- 
served ones given in Eq. (3) through an iterative Ncwton- 
Raphson method. Note that the component associated 
to Y e in the system (4) of hydrodynamic evolution equa- 
tions is treated as a passive advection equation, which 
does not contribute to the characteristic structure in the 



form of eigenvalues and eigenvectors required by some 
flux solvers. 

To numerically solve the metric equations we utilize 
an iterative nonlinear solver based on spectral methods. 
The spectral grid of the metric solver is split into 6 radial 
domains with 33 radial and 17 angular collocation points 
each. The combination of HRSC methods for the hy- 
drodynamics and spectral methods for the metric equa- 
tions (the Mariage des Maillages or "grid wedding" ap- 
proach) in a multidimensional numerical code has been 
presented in detail in [62]. Even when using spectral 
methods the calculation of the spacetime metric from 
the system (12-14) of elliptic equations is computation- 
ally expensive. Hence, in our simulations the metric is 
updated only once every 100/10/50 hydrodynamic time 
steps before/during/after core bounce, and extrapolated 
in between. The numerical adequacy of this procedure 
has been tested and discussed in detail in [14]. 

In this study we focus on the gravitational wave 
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signal associated with core bounce. As demonstrated 
by [17, 25], effects that may break rotational symmetry 
are most likely unimportant in this context. Hence, we 
assume axisymmetry and in addition impose symmetry 
with respect to the equatorial plane. 

The CoCoNuT code utilizes Eulerian spherical coor- 
dinates {r, 9}, and for the computational grid we choose 
250 logarithmically-spaced, centrally-condensed radial 
zones with a central resolution of 250 m and 45 equidis- 
tant angular zones covering 90°. A small part of the 
grid is covered by an artificial low-density atmosphere ex- 
tending beyond the core's outer boundary defined where 

p < io-Vc,i- 

We also note that we have performed extensive res- 
olution tests with different grid resolutions to ascertain 
that the grid setup specified above is appropriate for our 
simulations. 



IV. COLLAPSE DYNAMICS AND WAVEFORM 
MORPHOLOGY 

A. Generic waveform type 

We begin our discussion with an analysis of the gravi- 
tational radiation waveform emitted during core bounce 
as an indicator for the influence of the EoS, the progen- 
itor structure, and the precollapse rotational configura- 
tion on the collapse and bounce dynamics. In Fig. 2, we 
present example waveforms for representative collapsing 
cores selected from the investigated parameter space of 
models (i.e., less or more massive progenitors with slow 
or rapid precollapse rotation, varying degree of differen- 
tial rotation, and using either the Shen EoS or LS EoS). 
The waveforms of all models arc of type I, hence ex- 
hibit a positive prebounce rise and a large negative peak, 
followed by a ring-down. In the light of a considerably 
extended parameter space in terms of EoS and progeni- 
tor mass of the rotating core collapse models investigated 
in this work, we thus confirm the observation presented 
in [17, 18, 38, 41] that in general relativistic gravity 
all models with microphysics exhibit gravitational wave 
burst signals of type I. 

As already inferred in [18, 41], this signal type can be 
classified into three subtypes, which, however, do have in 
common the same qualitative features of a type I wave- 
form described above: 

(1) For a slowly rotating core, prompt convective over- 
turn at early postbounce times after the pressure- 
dominated bounce adds a significant low-frequency 
contribution to the regular ring-down signal (see, 
e.g., model sllAlO07). 

(2) In the case of moderately rapid rotation which 
still leads to a pressure-dominated bounce, this 
convection is effectively suppressed due to the 
growing influence of angular momentum gradi- 
ents [33, 68] and docs not strongly stand out in 
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FIG. 2: Time evolution of the gravitational wave amplitude 
h for representative models with different precollapse rota- 
tion profiles using the Shen EoS (red lines) or LS EoS (blue 
lines). Models with slow and almost uniform precollapse ro- 
tation (e.g., S11A1O07) develop considerable prompt post- 
bounce convection visible as a dominating lower-frequency 
contribution in the waveform, while the waveform for both 
models with moderate rotation (e.g., sllA3013, sl5A2O05, 
s20A2O09, S40A1O07, or el5a) and rapidly rotating models 
which undergo a centrifugal bounce (e.g., s40A3O13 or e20b) 
exhibit an essentially regular ring-down. Time is normalized 
to the time of bounce tb- 



the postbounce ring-down signal (see, e.g., mod- 
els sllA3013, S15A2O05, s20A2O09, S40A1O07, or 
el 5a). 

(3) If rotation is sufficiently rapid, the core bounces at 
subnuclear or only slightly supernuclear densities 
due to the increased effects of centrifugal forces, re- 
flected by a significant widening of the bounce peak 
of the waveform and an overall lower frequency of 
the signal (see, e.g., models s40A3O13 or e20b). 
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FIG. 3: Collapse dynamics of all investigated models in the 
parameter space of precollapse rotational configuration (spec- 
ified by the precollapse angular velocity S7 Ci i at the center and 
the precollapse differential rotation length scale A), progen- 
itor mass Mp r0 g, and EoS. Models marked by unfilled/filled 
circles undergo a pressure-dominated bounce with/without 
significant early postbounce convection, while models marked 
with crosses show a single centrifugal bounce. The EoS is 
encoded as in Fig. 1, while small/medium/large symbols rep- 
resent the precollapse rotation parameter A1/A2/A3. For 
better visibility, the symbols for the same M pr0 g but different 
EoS are spread a bit in the vertical direction. Note also that 
in this and the following plots that encode the precollapse 
rotational configuration in the form of the parameter fii, we 
refrain from including models el5a, el5b, e20a, and e20b, as 
these have a precollapse rotation profile that is not given by 
the analytic rotation law (18). 



Fig. 2 also demonstrates that for comparable precollapse 
rotational configuration (as specified by the parameters 
A and 17;) the impact of the EoS on the collapse dynam- 
ics and, hence, the gravitational wave signal is small. In 
Table 3 we mark each model with its type of collapse dy- 
namics, and in Fig. 3 we encode that type in the parame- 
ter space spanned by rotational configuration, progenitor 
mass/model, and EoS. 

For our set of collapse models, only in four cases (mod- 
els S11A1013, S15A1O07, S20A1O09, and s40A1005) the 
LS EoS yields a signal with dominant convective contri- 
bution while the Shen EoS docs not, and only a single 
model (sl5A2015) changes its collapse type from a cen- 
trifugal bounce to a pressure-dominated bounce when re- 
placing the Shen EoS with the LS EoS. However, Fig. 3 
shows that the transition between the three different col- 
lapse and waveform subtypes occurs at different precol- 
lapse rotational configurations for the various progenitor 
masses. This is a consequence of differences in the mass 
Micb of the inner core at bounce as discussed in Sec- 
tion IV C. 

The growth of the strong prompt early postbounce 
convection in slowly rotating models depends sensitively 



FIG. 4: Mass M- lc ^ of the inner core at the time of bounce for 
all models versus the precollapse initial central angular veloc- 
ity fi c .i. The progenitor model, the EoS, the initial rotation 
parameter A, and the collapse dynamics are encoded as in 
Fig. 3. 



on the seed perturbations resulting from the numer- 
ical scheme/grid. In nature, prompt convection will 
be influenced by random small-scale to large-scale vari- 
ations in the final stages of silicon burning that are 
frozen in during collapse. We point out that the dura- 
tion of the prompt postbounce convection is most likely 
overestimated in our approach, since in full postbounce 
radiation-hydrodynamics calculations, energy deposition 
by neutrinos in the immediate postshock region rapidly 
smoothes out the negative entropy gradient left behind 
by the shock (see, e.g., [20, 69]) and quickly damps this 
early convective instability. 



B. Influence of general relativity and 
deleptonization 

The general type of collapse and bounce dynamics of 
the core, i.e., pressure-dominated or centrifugal bounce, 
can be influenced (provided that the description of grav- 
ity and neutrino effects are identical) by the progeni- 
tor core stratification and thermodynamic structure, the 
amount and precollapse distribution of angular momen- 
tum, and the properties of the EoS in the density regime 
just below the stiffening threshold [18, 41]. These condi- 
tions influence the mass Mj Cj b of the homologously con- 
tracting inner core at bounce, which in turn determined 
the region that is dynamically relevant at bounce and 
sets the initial size of the proto-neutron star. 

In Fig. 4 we show the resulting variation of M; Cj b with 
O c ,i; progenitor model, precollapse differential rotation 
parameter A, EoS, and collapse type (encoded via sym- 
bols as in Fig. 3). The details of the variation of Mj Ci b 
with progenitor, EoS and rotational configuration will be 
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TABLE III: Summary of relevant quantities from the rotating collapse of all iron core models. p max ,b is the maximum density 
in the core at the time of bounce, |/i| max is the peak value of the gravitational wave amplitude for the burst signal, while /3b 
and /3 p b are the rotation rates at the time of bounce and late after core bounce, respectively. Models marked by unfilled/filled 
circles undergo a pressure-dominated bounce with/without significant early postbounce convection, while models marked with 
crosses show a single centrifugal bounce. The values left/right of the vertical separator (|) are for the Shen/LS EoS. 
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discussed in Section IV C. In the following, without loss 
of generality, we focus on a single progenitor and discuss 
the influence of general relativity and deleptonization on 
the collapse dynamics and the gravitational wave burst 
signal along the lines of the discussion in [18, 41]. 

In order to assess the individual influence of rclativis- 
tic effects and deleptonization, and to explain the ab- 
sence of type II and III burst signals in microphysical 
general relativistic models, in [18, 41] we compared col- 
lapse models of the s20 progenitor using the Shen EoS 
and a description for deleptonization with models using 
a simple hybrid polytropic/7-law EoS. We selected the 
adiabatic index 7 cos of these simple models in such a way 



that the transition between pressure-dominated bounce 
and centrifugal bounce occurs at the same precollapse 
rotation rate $ as for the microphysical models. With 
this method we were able to demonstrate that the influ- 
ence of deleptonization can be approximated by a cor- 
rection Ajv ~ 0.03 that must be applied to the estimate 
of the average EoS adiabatic index 7 e0 s,Shen — 1-32 in 
the density interval between 10 12 and 10 14 g cm' 3 . This 
leads to a generic value for the effective adiabatic index 
7cff,Shon = 7cos,shcn - — 1-32 - 0.03 = 1.29, prac- 
tically independent of the precollapse rotational configu- 
ration, both in Newtonian gravity and general relativity 
(where relativistic effects are accounted for by a correc- 
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FIG. 5: Boundary between pressure-dominated and centrifu- 
gal bounce in the j e s-/3i plane for s20 progenitor models using 
the hybrid EoS in Newtonian gravity (dashed lines) and gen- 
eral relativity (solid lines). The curved dotted lines show the 
Newtonian results shifted by — A7 gr = 0.015. The transi- 
tion points for models using the microphysical EoS without 
and with deleptonization, again for Newtonian gravity (cir- 
cles) and general relativity (bullets), lie in the shaded areas 
around 7eos,Shcn — 1-32 and 7 e fi,shen — 1-29, respectively, for 
the Shen EoS (top panel) and around 7 CO s,ls — 1.3225 and 
7efi,LS — 1.285, respectively, for the LS EoS (bottom panel). 



tion A7gr ~ —0.015). A graphic representation of this 
argument is shown in the top panel of Fig. 5, that is 
identical to Fig. 2 in [18] and Fig. 4 in [41], and which 
we include here for completeness. 

The estimate 7 a ~ 1.29 for microphysical models also 
allows us to explain the suppression of multiple centrifu- 
gal bounces with an associated type II waveform in a 
straightforward way, since this type of collapse occurs 
only in the respective hybrid EoS models with an ef- 
fective adiabatic index that is much closer to 4/3, i.e., 
7cff > 1.31. Rapid collapse dynamics that is character- 
ized by a type III burst signal is also not realized in mi- 
crophysical models of massive star collapse, as it requires 
a mass of the inner core at bounce Mi Cj b < 0.3 Mq [13] 
which is considerably smaller than those found in micro- 
physical models with any of our progenitors, for which 
we find Mi C > 0.4 Mq (see Fig. 4, and also the discussion 
in Section IV C). However, in [25, 70] is was suggested 
that rapid collapse dynamics and a type III burst signal 
may be associated with very efficient electron capture in 




FIG. 6: Adiabatic index 7 cos of the Shen EoS (red line) and 
LS EoS (blue line) versus the maximum density p max in the 
collapsing core for model s20A2O09. Although p max , which 
for this model is located in the center of the core, does not 
follow a trajectory of constant specific entropy, s is still ap- 
proximately conserved in the prebounce phase. Inset: Mag- 
nified view of 7 cos in the dynamically most relevant density 
range between 10 12 and 10 14 g cm -3 . The average value of 
7cos in this density regime is roughly 1.32 for both EoSs. 



the accretion-induced collapse of massive, rapidly rotat- 
ing white dwarfs. 

Finally, we point out that calculations with 7 er r used in 
the hybrid EoS have the tendency to underestimate the 
mass Mj Ci b at bounce compared to the full microphys- 
ical treatment. This is a consequence of the fact that 
in these calculations 7 e ff is kept constant throughout the 
collapse, leading to a reduction of the inner core mass 
Mj c already at much earlier collapse stages than in mi- 
crophysical models. The underestimated Afj Cj b, in turn, 
leads to gravitational wave burst signals from bounce 
in those simple models that are quantitatively or even 
qualitatively incorrect (as in the case of type III signals 
which do not occur in microphysical models). Hence, 
while useful for understanding the collapse dynamics, the 
7eft approach cannot replace the full microphysical treat- 
ment with a nonzero-temperature microphysical EoS and 
deleptonization as performed in the present work. 



C. Influence of the equation of state and 
progenitor model 

At densities below p nuc the total fluid pressure is dom- 
inated by the contribution from the degenerate electrons, 
hence the two microphysical EoSs should lead to rather 
similar dynamics in the infall phase of collapse. This is 
also reflected in the very similar behavior of their adia- 
batic indices 7 cos as shown in Fig. 6. 

In the bottom panel of Fig. 5 we demonstrate that the 
same influence of general relativistic effects and delep- 
tonization as discussed in the previous Section IV B ap- 
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FIG. 7: Boundary between pressure-dominated and centrifu- 
gal bounce in the 7 e ff— fi c ,i plane for models of all progenitors 
using the hybrid EoS in general relativity. The transition 
points for models using the microphysical EoS with delep- 
tonization (bullets) lie in the shaded area around 7„g ~ 1.28. 
Except for the rotation profile A2 of the sl5 progenitor, the 
locations of these points are identical for the two EoSs. Note 
that for the Al profiles of any progenitor (and the A2 profile 
for the sll progenitor) we do not observe a centrifugal bounce 
for any value of fi c ,i. In contrast to Fig. 5 we use here the 
fi c ,i instead of j3i as parameter for the precollapse rotational 
configuration (see the discussion in Section HE). 



plies for the s20 progenitor when the LS EoS is used 
instead of the Shen EoS. We obtain values of 7 0O s,ls — 
1.3225 for the adiabatic index of the EoS (without delep- 
tonization) and 7eff,LS — 1-285 for the effective adiabatic 
index (including deleptonization), which is in very close 
agreement with the values deduced from the simulations 
using the Shen EoS. As shown in Fig. 7, now only for 
general relativistic gravity, there is some spread of 7 e ff 
with progenitor mass/model, but on average, we find 
7 ff ~ 1.28 for the sll, sl5, and s40 progenitor models. 

Again following the line of arguments presented in [18, 
41], the combination of a low effective adiabatic index 
7ofi < 1.31 and a high inner core mass M- lc > 0.4 M & at 
bounce results in a type I gravitational wave burst signal 
for all our models, independent of the EoS or progen- 
itor model. Note that creating Figs. 5 and 7 we have 
performed additional simulations of microphysical mod- 
els that are more narrowly spaced in Q- hC (and corre- 
spondingly in /3{) than the ones listed in Table II. As 
a result, those figures reveal a small dependence of the 
transition between pressure-dominated bounce and cen- 
trifugal bounce (i.e., the location of the bullets and cir- 
cles in the direction of the abscissa) on the EoS, which is 
generally not apparent from Table II. 

Although the sensitivity of the deleptonization and 
collapse dynamics on the progenitor mass and EoS is 
only small, Fig. 4 still reveals a dependence of the in- 
ner core mass Mi Cj b at bounce both on the EoS and (in 
particular) on rotation. Furthermore, Mj Ci b varies non- 
monotonically with the progenitor mass M prog . In the 
absence of rotation, M; c .b is solely determined by the 



mean trapped lepton fraction Yj cp = Y e + Y v and specific 
entropy s in the inner core [1, 2, 27, 71] with a roughly 
quadratic dependence on both quantities. Since we em- 
ploy the same Y e (p) parametrization (based on model 
s20) for all models with the same EoS, the variations in 
Mic,h are caused by differences in the specific entropy 
in the precollapse iron core. Generally, the specific en- 
tropy in the iron core increases with progenitor mass, 
but, in particular in the mass range of ~ 18-25M Q , 
the relationship of progenitor mass and specific core en- 
tropy can be non-monotonic (see, e.g., [58]). However, 
note that the systematics for Mj Cj b with progenitor mass 
seen in Fig. 4 are possibly overemphasized by our Y e (p) 
parametrization and may be less pronounced in full ra- 
diation transport simulations which remain to be carried 
out in the future. 

For a rotating collapse, the variations of Mj Cj b with 
progenitor mass are amplified, while the overall system- 
atics are preserved. Obviously, a more massive and hence 
more extended inner core is more susceptible to the in- 
fluence of centrifugal forces (which scale proportional to 
the radius r) than a less massive and thus more compact 
inner core. This behavior is confirmed by Fig. 4, which 
depicts the dependence of the mass Mi Cj b of the inner 
core at bounce on the precollapse central angular veloc- 
ity fi Cj i, the EoS, and the differential rotation parameter 
A. Models with comparably large precollapse iron core 
specific entropy (and also large iron core mass) and, thus, 
larger Mj Cj b already in the nonrotating case, show a more 
pronounced increase of Mj c .b with rotation than models 
with lower precollapse specific entropy (and also smaller 
iron core mass). The scaling of M- lc ^ with f2 c ,i, at fixed 
differential rotation parameter A, is linear for small to 
intermediate fl Cj i and turns approximately quadratic for 
the most rapidly rotating configurations. On the other 
hand, when increasing the degree of differential rotation 
A at fixed f2 Cl i, Mj Ci b decreases since then the angular 
velocity in the outer parts of the inner core and conse- 
quently centrifugal support drops. 

We also observe that the impact of the EoS on the 
mass of the inner core manifests itself only via an al- 
most constant positive relative increase in Mj Cj b of ~ 10% 
when changing from the LS EoS to the Shen EoS, prac- 
tically independent of rotation and progenitor mass (see 
Fig. 4). Again, the mean electron (respectively lepton) 
fraction and specific entropy in the inner core are the key 
to understanding these systematics. The representative 
s20 progenitor model used to parametrize Y e {p) in this 
study yields minima for Y e in the center of the core at 
bounce of ~ 0.249 and ~ 0.241 for the Shen EoS and the 
LS EoS, respectively. This relative difference of ~ 3.3% 
translates into a difference in Mj Ci b of ~ 7%, assuming 
that the mass of the inner core scales quadratically with 
Y e , which slightly underestimates the actual change. We 
attribute the remaining difference to variations in the 
specific entropy s of the inner core at bounce due to the 
slightly more efficient electron capture in the models with 
the LS EoS. 
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FIG. 8: Maximum density p maXj b in the star at the time of 
bounce for all models versus the precollapse central angular 
velocity fl c ,i- For moderately fast or rapidly rotating models, 
which do not exceed nuclear density at bounce, p max ,b is al- 
most identical for the Shen EoS (dark hues) and the LS EoS 
(light hues), while for slowly rotating models the difference 
reaches Ap max> b — 10 14 g cm -3 in the nonrotating limit. The 
progenitor mass, the EoS, the precollapse rotation parameter 
A, and the collapse dynamics are encoded as in Fig. 3. 



We point out that the progenitor models el5a, el5b, 
c20a, and e20b, which already come with a rotation pro- 
file from the stellar evolution calculation, are very well 
represented in terms of collapse dynamics, waveform, 
and postbounce rotation state by members of our model 
set with an artificially added precollapse rotation profile, 
specifically the models sl5A2O09, sl5A2015, s20A2O09, 
and s20A2O15. For this reason we refrain from separately 
discussing those special models in the entire paper. 



D. Influence of differential rotation 

Increasing the degree of differential rotation by lower- 
ing the value of the differential rotation parameter A at 
fixed precollapse central angular velocity r2 c ,; results in 
less centrifugal support in outer core regions and, as al- 
ready pointed out in Section IV C, in a smaller mass Mj Cj b 
of the inner core at bounce. Consequently, a higher value 
of Cj i is necessary for a stronger differentially rotating 
precollapse core to become significantly affected by cen- 
trifugal forces during the collapse. This is confirmed by 
Fig. 7 which displays the systematics of the transition be- 
tween pressure-dominated and centrifugal bounce for our 
set of progenitors and the A2 and A3 rotation profiles. 
Compared to the transition values of f2 C; i for the A2 pro- 
file, the A3 profile requires a roughly 20-40% higher f2 Cji 
(varying slightly with progenitor model) for a transition 
from pressure-dominated to centrifugal bounce. 

In previous extensive parameter studies of rotating 
stellar core collapse (see, e.g., [14-16]) the effect of dif- 



FIG. 9: Relative change Ap maXj b,rci of the maximum density 
at bounce when changing from the Shen EoS to the LS EoS 
for all models versus the precollapse central angular velocity 
Q c ,i- The progenitor mass and the precollapse differential ro- 
tation parameter A are encoded as in Fig. 3, while the collapse 
dynamics are not specified. 



ferential rotation was studied in model sequences in the 
parameter space spanned by the precollapse differential 
rotation parameter A and the precollapse rotation rate 
/?;. At a constant /3i, more differentially rotating models 
require a larger f2 Cj ; than less differentially rotating ones 
and experience core bounce at lower densities. Hence, 
at fixed /3j, more differentially rotating models are gener- 
ally more affected by centrifugal effects. Our s20 model 
series is constructed as a sequence of fixed # for each 
of the rotation profiles Al, A2, and A3 (see Table II), 
and therefore permits a direct comparison with preceed- 
ing work. Our results confirm qualitatively the previ- 
ously identified systematics (see Table III). However, in 
contrast to more simplistic simulations, the combination 
of general relativity and deleptonization in our models 
weakens the overall impact of centrifugal effects on the 
collapse dynamics (see Section IV B), and consequently 
leads to much smaller quantitative changes in the charac- 
teristic collapse variables (such as p m ax,b) |^|max, -^ic.b, 
or /3b) when varying the degree of differential rotation. 



V. STRUCTURE OF THE POSTBOUNCE CORE 
AND IMPACT ON THE WAVE SIGNAL 

A. Equation of state at supernuclear densities and 
maximum density in the core 

From Table III it is apparent that the change from 
the Shen EoS to the LS EoS in an otherwise identical 
model results systematically in an increase of the peak 
maximum density p maXib at bounce, i.e., /9 max ,b,LS > 
Pmax,b,Shen- This result is in agreement with the previ- 
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ous Newtonian study of Kotake et al. [32] who compared 
simulations of a single model carried out with the Shen 
EoS and the LS EoS. 

For centrifugally bouncing models, which only 
marginally exceed or even remain below p n uc at bounce, 
the absolute increase in the maximum core density at 
bounce is small, exhibiting a maximum Ap max .b = 

Pmax,b,LS ~ Pmax,b,Shcn = 0.13 X 10 14 g Cm" 3 for 

model s40A3O09 (leaving aside the exceptional models 
s40A2O13 and s40A2O15 which we will discuss sepa- 
rately later). This is another manifestation of the similar- 
ity of the two microphysical EoSs at subnuclear densities 
(see also Section IV C). 

For slowly or at most moderately fast rotating mod- 
els that undergo pressure-dominated bounce and whose 
center exceeds supernuclear density at (and also after) 
bounce, Ap max b can amount up to 1.19 x 10 14 g cm -3 
for model sllAlOOl, the most slowly rotating model of 
the sll model scries. This strong impact of the EoS 
can be readily explained by the fact that at supernuclear 
densities the LS EoS is considerably softer than the Shen 
EoS. Fig. 6 shows a difference in the adiabatic index 7 eos 
between the two microphysical EoSs (for a representa- 
tive model) of about A7 eos ~ —0.5 at those densities, 
where nuclear forces have an essential impact on the EoS 
properties. 

The large effect of the EoS seen in p max ,b in models 
where this quantity exceeds p nuc does not contradict our 
observation that the EoS has little impact on the col- 
lapse dynamics, since once the core plunges into the su- 
pernuclear density regime, where stronger differences in 
the two microphysical EoSs emerge, the mass Mj c .b of 
the inner core at bounce is already fixed and the bounce 
dynamics (pressure-dominated or centrifugal) is already 
determined. 

The impact of the EoS on p m ax.b is also visualized in 
Fig. 8. As expected, for moderately or rapidly rotating 
models, whose central parts do not reach high supernu- 
clear densities at bounce, the difference in p max ,b grad- 
ually decreases. Fig. 8 also reveals that the two models 
s40A2O13 and s40A2O15 (marked by two dark and light 
blue crosses at intermediate values of fi c ,i, respectively) 
are the ones that undergo a clear centrifugal bounce for 
both EoSs with the lowest value of fl c .i of all models with 
the A2 rotation profile. 

The convergence of p max ,b for the two microphysical 
EoS with increasing rotation can also be observed in the 
relative difference Ap maXib , r ci = p ra ax,b,Ls/Pmax,b,shcn - 1 
shown in Fig. 9, which starting from Ap max .b,rci — 
35-40% in the nonrotating limit first declines linearly 
with f2 Cj i until it levels off at roughly constant val- 
ues. However, the largest values are obtained with 
Ap maX) b,roi — 40% and 48% for the rapidly rotating and 
centrifugally bouncing models s40A2O13 and s40A2O15, 
emphasizing their exceptional nature. This particular 
behavior results from a combination of two effects, ex- 
hibited by only these two models in our entire model 
set. First, when switching from the Shen EoS to the LS 
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FIG. 10: Time evolution of the maximum density p max for 
representative models with different precollapse rotation pro- 
files using the Shen EoS (red lines) or LS EoS (blue lines). 
While models with at most moderate precollapse rotation 
rates (e.g., sllAlO07 or s20A2O09) undergo a pressure- 
dominated bounce at supernuclear densities, rapidly rotating 
models (e.g., e20b) experience a centrifugal bounce. 



EoS the inner core mass Mi c> b at bounce significantly de- 
creases (see Fig. 4). Hence in the LS EoS variant the two 
models experience weaker rotational support (in partic- 
ular with the differential rotation parameter A2; see also 
Section IV C). Second, the two models bounce in a den- 
sity regime (see Table III and Fig. 6) where the LS EoS 
exhibits a smaller 7 oos than the Shen EoS, resulting in 
less pressure support when the LS EoS is used. The 
combination of weaker rotational support and pressure 
support when using the LS EoS can then readily explain 
the excess in p max ,b,LS compared to p max ,b,shcn in the two 
exceptional s40 models. 

A higher value of the maximum density p ma x in the 
collapsed core for the LS EoS is not limited to the time 
of bounce, but typically also remains in the nascent 
proto-neutron star at later postbounce times, as shown 
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FIG. 11: Peak value |/i| max of the gravitational wave ampli- 
tude at 10 kpc distance for the burst signal (neglecting possi- 
bly larger contributions from postbounce convection at later 
times) for all models versus initial precollapse central angular 
velocity £l c ,i- The progenitor mass, the EoS, the precollapse 
differential rotation parameter A, and the collapse dynamics 
are encoded as in Fig. 3. 



FIG. 12: Relative change A|fe| max , ro i of the peak value |/i| ma x 
of the gravitational wave amplitude for the burst signal when 
changing from the Shen EoS to the LS EoS for all models ver- 
sus the precollapse central angular velocity Q c ,i- The progen- 
itor mass and the precollapse differential rotation parameter 
A are encoded as in Fig. 3, while the collapse dynamics are 
not specified, as for some models it depends on the EoS. 



in Fig. 10 for models representing the three collapse type 
and waveform subclasses (see Section IV A). Only very 
rapidly and thus centrifugally bouncing models such as 
model e20b in Fig. 10 have a time evolution of /3 max that 
is practically independent of the chosen EoS. 

We point out that in our discussion we always make use 
of the maximum density p max instead of the central den- 
sity p c , since, after bounce, some of the most rapidly ro- 
tating and thus centrifugally bouncing models develop a 
slightly toroidal density structure with an off-center den- 
sity maximum that is at most 20% larger than p c . This 
is much less extreme than for models with the simplified 
hybrid EoS treatment, where the maximum density was 
found to be several orders of magnitude larger than the 
central density in extreme cases [13, 14]. 



B. Structure of the postcollapse core and peak 
waveform amplitude 

Since the LS EoS leads to higher central densities at 
bounce, one can on the one hand expect higher gravita- 
tional wave peak amplitudes in the burst signal from core 
bounce, as a denser and more compact core should yield 
in an increase of the contribution to the quadrupole mo- 
ment from the central parts of the core. Furthermore, the 
associated shorter dynamical times also lead to a more 
rapid time variation in the quadrupole formula. On the 
other hand, the higher compactness of the inner core of 
a model run with the LS EoS results in lower densities at 
intermediate and large radii than in the less compact core 
of the corresponding model with the Shen EoS. In turn, 



this may lead to an effectively smaller total quadrupole 
moment and thus to a decrease in the signal amplitude 
compared to the counterpart model with the Shen EoS. 
We now test which of these two competing effects domi- 
nates in our models. 

In Fig. 11 we show the peak value |/i| m ax of the grav- 
itational wave amplitude for the burst signal from core 
bounce (see also Table III), where we neglect any possibly 
larger contributions at later times for models with strong 
prompt postbounce convection. For slowly or at most 
moderately rapidly rotating cores, |/i| max rises steeply 
with increasing f2 Ci j and covers a range of more than two 
orders of magnitude for our selection of initial models. 
For rapid rotation, when centrifugal forces become dy- 
namically important and can be the dominant factor at 
bounce, the peak amplitude |/i| max saturates and even 
decreases again at very high O c ,i. This behavior is a 
consequence of centrifugal support, which prevents such 
rapidly spinning cores from reaching high densities and 
more extreme compactness as well as being subject to 
short variations of the quadrupole moment (see also the 
discussion in Section VII and in [16]). 

For each precollapse rotational configuration (i.e., at 
constant O c i and differential rotation parameter A in 
Fig. 11), the value of |/i| max depends indirectly on the 
mass Mp rog of the progenitor via the mass Mj Ci b of the in- 
ner core at bounce. As already discussed in Section IV C, 
A/ic.b does not depend in a monotonic way on Af prog , 
but for our standard model set increases in the order of 
the progenitor models sll, s20, sl5, and s40. Therefore, 
for pressure-dominated bounce models the amplitude of 
the gravitational wave signal directly scales with M- lc ^ in 
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FIG. 13: Radial profiles of the density p c in the equatorial 
plane at the time of bounce for model s20A3O09 (top panel) 
and model s40A2O13 (bottom panel) using the Shen EoS (red 
line) and LS EoS (blue line). In the central parts of the proto- 
neutron star (for these models at radii smaller than the cross- 
ing radius r ~ 5.5 km and r ~ 23.5 km, respectively, marked 
by dotted lines) the LS EoS leads to higher densities. 

the obvious sense that more massive inner cores produce 
stronger gravitational wave emission. 

What cannot be extracted from Fig. 11 is a clear 
effect of the choice of the EoS on |/i| max , despite the 
strong difference in /o max ,b we observe between models 
using the Shen EoS and the LS EoS. When plotting 
the relative change A|/i| maXjre i = |/i| m ax,Ls/Nmax,shon - 
1 obtained by changing from the Shen EoS to the 
LS EoS for the same initial model (as presented in 
Fig. 12), the majority of models shows a decrease of 
Hmax,LS compared to |/i| max ,shcn- Only six out of the 
68 models (sllA2O05, sllA2O07, sllA3O05, S20A1O05, 
S20A1O07, and s40A3O09) listed in Table III exhibit a 
larger |/i| max when the LS EoS is used. 

This behavior is similar to the situation discussed by 
Dimmelmeier et al. [14] who compared collapse dynamics 
and gravitational waveforms obtained from Newtonian 
and general relativistic collapse simulations with the sim- 
ple hybrid EoS. They showed that for |/i| max the global 
density distribution in the core at bounce is decisive, not 
the local maximum density value. In their simulations, 
the general relativistic variants consistently produced an 
increase of /0 max .b compared to their Newtonian counter- 
parts. Still, they found that the peak value |/i| max of 
the gravitational wave amplitude actually decreases for 
most models when general relativistic effects are taken 
into account. 

In [14], the negative A|/i| maXjre i observed in many mod- 
els when comparing Newtonian and relativistic simula- 
tions could be attributed to the "density crossing" that 
occurs at some radius inside the inner core at bounce: 




r [km] 

FIG. 14: Radial profiles of the weighted density p e r 2 in the 
equatorial plane at the time of bounce for model s20A3O09 
(top panel) and model s40A2O13 (bottom panel) using the 
Shen EoS (red lines) and LS EoS (blue lines). The vertical 
lines mark the crossing radius. 

The general relativistic simulation of a model yields a 
higher density inside that (angle-dependent, due to ro- 
tation) radius, while for larger distances from the cen- 
ter, p is smaller compared to the Newtonian simulation. 
Here, we vary the EoS rather than the description of 
gravity, but we observe a very similar density crossing in 
models that show |/i| max ,LS < Mmax.Shen' In Fig. 13 we 
demonstrate this for models s20A3O09 (representative 
for a pressure-dominated bounce) and s40A2O13 (repre- 
sentative for a centrifugal bounce). 

Following the argument in [14], we plot the weighted 
density pr 2 in Fig. 14, since this is the relevant quantity 
in the integrand of the quadrupolc gravitational wave 
formula. Although the larger pr 2 of the model with the 
LS EoS gives a larger quadrupolc contribution out to the 
crossing radius, in most models the larger pr 2 in the outer 
parts of the core in the variant with the Shen EoS more 
than compensates this and ultimately leads to a larger in- 
tegral quadrupolc moment and, thus, to a stronger grav- 
itational wave burst. We note that in [14], all models 
whose collapse type did not change exhibited lower peak 
waveform amplitudes (A|ft,| maX)re ] < 0) when going from 
Newtonian to general relativistic gravity. In contrast, go- 
ing from the relatively stiff Shen EoS to the softer LS EoS 
results in a less clear trend with a few models exhibiting 
A|/i| max ,rci > 0. This suggests a less dramatic impact of 
a change from the Shen EoS to the LS EoS compared to 
altering the description of gravity. 

For the small set of our models with A|/i| maXire i > 
we are neither able to identify any obvious and system- 
atic correlation with model parameters nor do we find 
any clear systematics of A|/i| max rc i > with Ap max b, rel- 
it appears that the sign and magnitude of A|/i| maXire i 
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depends sensitively and in a complicated way on the de- 
tails of the collapse dynamics in each individual model. 
Hence, we can explain only why for a specific model dif- 
ferences in the density structure at bounce between the 
model variants with the Shen EoS and the LS EoS lead 
to an observed A|/i| maXire i, but cannot predict A|/i| max re i 
based on precollapse model parameters. 



C. Frequency spectrum of the waveform and 
variation with the equation of state 

In contrast to the somewhat ambiguous impact of the 
EoS on the peak waveform amplitude, the effect of re- 
placing the Shen EoS with the LS EoS on the waveform 
peak frequency is unequivocal for models undergoing a 
pressure-dominated bounce. The increase in the maxi- 
mum density at bounce in the models with the LS EoS 
always results in a shift of the main peak in the waveform 
spectrum to higher frequencies. In the center panel of 
Fig. 15, we plot the waveform spectrum (i.e., the Fourier 
transform h of h) for model s20A3O09 as a representa- 
tive pressure-dominated bounce model. The spectrum of 
this model exhibits a distinct and narrow high-frequency 
peak at / m ax,Shcn = 710 Hz when using the Shen EoS, 
while the calculation of the same model with the LS EoS 
results in a peak at / m ax,LS = 744 Hz. Thus, for this 
particular model, the change in EoS shifts the frequency 
associated with the bounce peak by A/ max = +34 Hz. 
We observe similar values for A/ max in all models under- 
going pressure-dominated bounce. 

At frequencies below about 200 Hz, the waveform spec- 
trum of s20A3O09 exhibits a plateau which is due to 
the low-frequency contribution from prompt large-scale 
postbounce convection. Such a contribution is present 
in many models with slow to moderate rotation, but 
gradually decreases in magnitude and relevance with in- 
creasing rotation. As pointed out in Section IV A, our 
present numerical scheme has the tendency to overes- 
timate prompt postbounce convection compared to full 
radiation- hydrodynamics calculations. 

The waveform of the slowly spinning model s20AlO05, 
whose spectrum is shown in the top panel of Fig. 15, is 
dominated by such prompt postbounce convective mo- 
tions. Accordingly, for this model, there is a strong con- 
tribution to the spectrum at low frequencies, even ex- 
ceeding the still clearly discernible bounce peak at high 
frequencies. Nevertheless, also in this case the shift of 
the high-frequency bounce peak when replacing the Shen 
EoS by the LS EoS is obvious and obeys the systematics 
discussed above. 

With increasing rotation, centrifugal forces become 
more relevant and slow down the late phase of collapse 
and bounce. As a consequence, / max always retreats to 
lower frequencies. This is apparent in the spectrum of the 
centrifugal bounce model s20A3O13 shown in the bot- 
tom panel of Fig. 15. For this model, one can still iden- 
tify the high-frequency bounce peak, but now at signifi- 
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FIG. 15: Spectrum of the gravitational radiation waveform 
for model S20A1O05 (top panel), model s20A2O09 (center 
panel), and model s20A3O13 (bottom panel) using the Shen 
EoS (red line) and LS EoS (blue line), h is the Fourier trans- 
form in frequency space of the waveform amplitude h. The 
dotted lines mark the frequency / max at the maximum of the 
waveform spectrum, neglecting low-frequency contributions. 
The scale of the vertical axis is one order of magnitude per 
major tick mark. 

cantly lower frequencies around 400-500 Hz. Note that 
the low-frequency quasi-continuous part of the spectrum 
in centrifugally bouncing models such as s20A3O13 is due 
to rotationally slowed dynamics and stronger postbounce 
oscillations, and should not be confused with the low- 
frequency contribution from prompt convection in slowly 
rotating models. 

In Fig. 16, we plot / max for all models that undergo 
pressure-dominated bounce and thus exhibit a clearly 
visible high-frequency peak in their spectra that can be 
associated with the gravitational wave burst from core 
bounce. For all models the systematic increase of / max 
when changing from the Shen EoS to the LS EoS is appar- 
ent, and only for very few rapidly rotating models close 
to the threshold to centrifugal bounce the change of / max 
becomes small. In Table IV we summarize the arithmetic 
mean / max along with the respective absolute and relative 
differences between models using the Shen EoS and the 
LS EoS. Note that when computing / max we neglect the 
contribution below a cut-off frequency / cut = 250 Hz in 
order to exclude any influence from the possibly unphys- 
ically strong and prolonged early postbounce convection. 

In previous work [18], Dimmclmeier et al. discussed 
the detection prospects for the gravitational wave burst 
emitted in rotating core collapse models based on the s20 
progenitor and using the Shen EoS. To this end, they 
simulated a large set of models with varying precollapse 
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FIG. 16: Frequency / max at the maximum of the waveform 
spectrum for all models with a given progenitor mass versus 
the precollapse central angular velocity Q C) i. Only models 
which undergo pressure-dominated bounce are shown. The 
dotted lines mark the average / max when using the Shen EoS 
or LS EoS. The progenitor mass, the EoS, the precollapse 
differential rotation parameter A, and the collapse dynamics 
are encoded as in Fig. 3. 



rotation rates $ in the range from 0.05% to 4%, ap- 
proximately logarithmically spaced in 18 steps for each 
of the three rotation profiles Al, A2, and A3. For the 
current work, we have repeated the calculations of this 
model set (which is extended in terms of precollapse ro- 
tation compared to our standard models stated in Ta- 
ble II, but limited to one progenitor), this time with 
the LS EoS. While the models with the Shen EoS that 
undergo pressure-dominated bounce have an arithmetic 
mean peak frequency /max.shon ~ 718 Hz [18], we find 
/maxis ~ 758 Hz when using the the LS EoS. Thus for 
this particular model set the average relative frequency 
shift amounts to A/ maxre i ~ 5.6%. Both the average 
peak frequencies and their change with EoS are consis- 
tent with what we find for our standard model set using 
the four different progenitors and a more restricted vari- 
ety of precollapse rotation rates. 



TABLE IV: Average / max of the frequency at the maximum 
of the waveform spectrum for all models with a given progen- 
itor mass. A/ max and A/ maX)re i are the absolute and relative 
change of the frequency average, respectively, when changing 
from the Shen EoS to the LS EoS. 



Collapse 


/max,Shen 


/"max , L S 




^/max,rel 


model set 


[Hz] 


[Hz] 


[Hz] 


[%] 


sll 


733 


777 


44 


6.0 


sl5 


658 


702 


44 


6.7 


s20 


690 


724 


34 


4.9 


s40 


685 


716 


31 


4.5 



VI. DETECTION PROSPECTS FOR THE 
GRAVITATIONAL WAVE BURST SIGNAL 



In order to assess the detectability of the burst signal 
from core bounce, we compute the (detector-dependent) 
frequency-integrated characteristic signal frequency f c 
and dimcnsionless characteristic gravitational wave am- 
plitude h c using Eq. (22) and Eq. (23), respectively. We 
again exclude frequencies below 250 Hz from the inte- 
grals in an attempt to filter out dominant contributions 
from prompt postbounce convection in slowly rotating 
models. In Fig. 17, we plot h c against f c for the current 
LIGO detector [72] at a distance of 10 kpc. For compari- 
son with the detector sensitivity, we include its rms strain 
sensitivity curve. Note that the total energy emitted in 



gravitational waves ranges from E g 



3.5 x 10 



-10 



to 



5.3 x 10 8 in units of M Q c 2 (including the contribution 
from convection) for our standard models. 

The distribution of our standard set of models (as 
listed in Table II) in the h c -f c plane of Fig. 17 obeys 
straightforward systcmatics. The clustering in frequency 
of the large number of models undergoing a pressure- 
dominated bounce (marked by circles in Fig. 17) is obvi- 
ous. Very slowly rotating models, whose waveforms arc 
dominated by the imprint of prompt postbounce convec- 
tion (unfilled circles), exhibit the lowest values for h c , 
which increases with faster rotation (along arrow 1), re- 
flecting that the inner core at bounce becomes more mas- 
sive (cf. Sections VB and IV C). Despite the frequency 
cut at 250 Hz in the integral for f c , the low- frequency 
contribution from convection in the spectrum leads to an 
f c that is lower than the value obtained for more rapidly 
rotating models without significant postbounce convec- 
tion (filled circles). For the latter model class, h c simply 
grows with increasing precollapse rotation (along arrow 
2), now at practically constant f c . Even for these mod- 
els, f c is always lower than the average peak frequency 
/max of their waveform spectra, which amounts to 715 Hz 
for the 108 models of our standard model set (including 
the el5/e20 models) which exhibit a pressure-dominated 
bounce. This is a consequence of the detector character- 
istics, whose maximum sensitivity is at much lower fre- 
quencies between 100 and 200 Hz and thus accordingly 
lowers f c in comparison with a fiducial flat sensitivity 
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FIG. 17: Location of the gravitational wave burst signals from 
core bounce for all models (including the el5/e20 models) in 
the h c -fc plane relative to the sensitivity curves of the LIGO, 
assuming at a distance of 10 kpc. The meaning of the arrows 
1, 2, and 3 as well as area 4 are explained in the main text. 
The dotted line marks the average / max of the frequency at the 
maximum of the waveform spectrum. The progenitor model, 
the EoS, the initial rotation parameter A, and the collapse 
dynamics are encoded as in Fig. 3. 



curve. 

For rapid rotation, the influence of centrifugal forces 
on the collapse dynamics manifests itself as a centrifugal 
barrier that limits the characteristic amplitude h c (see 
also the discussion in Section VII A and Fig. 11). Simul- 
taneously, the characteristic frequency f c moves to in- 
creasingly lower values as faster rotation slows down the 
collapse (along arrow 3). Models that rotate so rapidly 
that they undergo a purely centrifugal bounce (marked 
by cross symbols in Fig. 17) constitute a practically sepa- 
rate class (area 4) in the h c —f c diagram somewhat below 
the maximum value of the amplitude h c , but at consid- 
erably lower frequencies f c . 

For very rapidly rotating models the imprint of cen- 
trifugal effects on various waveform characteristics (such 
as / m ax, fc, Hmax, or h c ) is quite pronounced and per- 
mits one to infer on the precollapse rotational config- 
uration in the case of a successful detection of gravita- 
tional waves from a core collapse event. As already noted 
in [18], in the case of moderate or slow rotation, which 
is the astrophysically most probable case [31, 59], the in- 
sensitivity of the waveform's frequency characteristics to 
variations in the precollapse configuration significantly 
obstructs the "inversion problem" of gravitational wave 
detection, i.e., the constraining of physical parameters of 
the precollapse core or of the nascent proto-neutron star 
from a detected waveform, leaving only the (e.g., maxi- 
mum or integrated characteristic) amplitude as an indi- 
cator of the rotational configuration. In addition, Fig. 16 
also implies that it will be very hard, if not impossible, 
to constrain other possibly unknown model parameters 



FIG. 18: Location of the gravitational wave burst signals 
from core bounce in the h c -f c plane relative to the sensitivity 
curves of various interferometer detectors (as color-coded) for 
an extended set of models with the progenitor s20 using the 
Shen EoS (dark hues) or LS EoS (light hues). The sources 
are at a distance of 10 kpc for LIGO, 0.8 Mpc for Advanced 
LIGO, and 15 Mpc for EURO. The dotted lines mark the av- 
erage /max of the frequency at the maximum of the waveform 
spectrum for the models when using the Shen EoS or LS EoS. 
Only the EoS and the collapse dynamics are encoded as in 
Fig. 3, but not the precollapse differential rotation parameter 
A. 



aside from rotation (such as EoS or progenitor mass) from 
the gravitational waveform of the burst signal from core 
bounce alone, since their effect on the burst waveform is 
small and no clear trends or systematics are discernible, 
which adds to the degeneracy of the inversion problem. 

As an example, we again single out the impact of the 
EoS on the waveform frequency while keeping the progen- 
itor model s20 fixed. For the particular, extended set of 
models with many precollapse rotation rates already dis- 
cussed in Section VC, we show in Fig. 18 the location of 
the waveform signals in the h c —f c plane for initial LIGO 
at a distance of 10 kpc, Advanced LIGO in broadband 
tuning [72] at a distance of 0.8 Mpc, and the projected 
EURO detector in xylophone mode [73] at a distance of 
15 Mpc (cf. Fig. 4 in [18]). All 54 s20 models of [18] 
using the Shen EoS along with the newly computed cor- 
responding models with the LS EoS are shown. 

It is obvious that the spread within the group of models 
with either the Shen EoS or the LS EoS is larger than the 
variation due to a change in the EoS, since the effect of 
the EoS on the characteristic signal frequency f c is small 
(comparable to A/ maXjro i, corresponding to a change of 
a few percent). The two EoSs considered here bracket 
the range from rather soft (LS EoS) to rather stiff (Shen 
EoS), and therefore it is unlikely that employing a larger 
variety of nonzero-temperature nuclear EoSs would lead 
to any more optimistic conclusions. 

Based on the relative positions of the models with re- 
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spect to the individual detector sensitivities, from Fig. 18 
we conclude (in agreement with previous work [17, 18]) 
that initial-LIGO-class detectors are sensitive only to 
signals coming from an event in the Milky Way, while 
Advanced-LIGO-class observatories could marginally de- 
tect events from other galaxies in the Local Group (e.g., 
M31 Andromeda at ~ 0.8 Mpc distance). For the pro- 
posed EURO detector in xylophone mode, we expect a 
very high signal-to-noise ratio (h c divided by the detector 
sensitivity at / c ). This detector could also observe many 
of the computed signals at a distance of 15 Mpc, i.e., in 
the Virgo cluster, for which one expects a favorably high 
event rate. 



VII. ROTATION OF THE PROTO-NEUTRON 
STAR 

The calculations presented in this study impose ax- 
isymmctry, hence are unable to track the development 
of rotationally-induced nonaxisymmetric structures and 
dynamics. Nevertheless, we can utilize the results from 
our simulations to assess the possibility of rotational 
triaxial instabilities during the collapse and early post- 
bounce phase. In this way we can (i) test the reliability 
of our present restriction to axisymmctry and (ii) put 
constraints on the relevance of the various types of such 
instabilities in a core-collapse event. 

Nonaxisymmetric rotational instabilities in proto- 
neutron stars have long been proposed as strong and po- 
tentially long-lasting sources of gravitational waves. In 
principle, the gravitational wave emission by a nonax- 
isymmetrically deformed proto-neutron star after bounce 
could easily exceed (see, e.g., [17, 25]) in total emitted 
energy (and, hence, in characteristic strain h c ) the gravi- 
tational wave burst from core bounce on which this paper 
is focussed. 

In the context of classical Newtonian theory of fluid 
equilibria (see, e.g., [74]), MacLaurin spheroids (i.e., ax- 
isymmctric, rigidly rotating, equilibrium configurations 
of uniform density) become unstable to nonaxisymmet- 
ric deformation when a nonaxisymmetric configuration 
with lower total energy exists at a given rotation rate /?. 
MacLaurin spheroids become dynamically unstable to de- 
formation into Riemann ellipsoids at (3 > f3d yn = 27%. 
At j3 > (3 SCC = 14% they become secularly unstable to tri- 
axial ellipsoidal deformation in the presence of dissipative 
processes (Jacobi ellipsoids via gravitational wave back- 
reaction known as the Chandrasekhar-Friedman-Schutz 
(CFS) instability [75, 76], or Dedekind ellipsoids via vis- 
cous processes). In both the dynamical and the secular 
case, the lowest-order deformation in terms of azimuthal 
nonaxisymmetric modes proportional to cxjp(imtp) is the 
to = 2 Kelvin (bar-) /-mode, where if is the azimuthal 
angle and the mode order to is an integer. 

Although Newtonian MacLaurin spheroids are highly 
idealized configurations, numerical studies (see, e.g., [77] 
and references therein) have shown that the above insta- 
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FIG. 19: Time evolution of the rotation rate (3 around the 
time of core bounce for various models of the s20 progenitor 
series computed with the Shen EoS at fixed precollapse degree 
A of differential rotation and varying the precollapse central 
angular velocity fl c ,i- Note that we have augmented this se- 
quence by three extra models s20A3O16 to s20A3O18 (with 
A = 3.00, 3.50, and 4.00, respectively) not listed in Table III. 

bility threshold /3dy n for the dynamical instability holds 
approximately even when differentially rotating com- 
pressible fluid configurations in general relativity are con- 
sidered. The situation may be different for the gravita- 
tional radiation back-reaction driven secular instability, 
since perturbative studies (see, e.g., [78]) predict an on- 
set at significantly lower (3 in general relativity than in 
the Newtonian case. However, fully relativistic nonlinear 
hydrodynamic studies of the secular instability remain 
yet to be carried out. 

Recently, a new kind of dynamical rotational nonax- 
isymmetric instability at a value of f3 much lower than the 
classical threshold has been discovered both in numeri- 
cal and perturbative studies (see, e.g., [17, 22, 25, 79-85] 
and references therein). This low-/3 instability (making 
the classical MacLaurin instability a "high"-/3 instabil- 
ity) appears to amplify nonaxisymmetric modes at points 
where their pattern speed a m (the eigenfrequency uj m di- 
vided by the azimuthal mode order to) coincides with the 
local angular velocity of the fluid [22, 81-83]. 



A. The rotational barrier in core collapse 

From first principles one can derive that the conser- 
vation of angular momentum during the collapse phase 
results in an increase of the angular velocity f2 of a 
representative Lagrangian mass element proportional to 
zu , where vj = r sin (9 is the distance from the rota- 
tion axis. Setting for simplicity w equal to the spheri- 
cal radial coordinate r (which, of course, only holds in 
the equatorial plane) , this translates into a scaling of the 
centrifugal force proportional to r -3 . The gravitational 
force, on the other hand, increases only like r~ 2 . Hence, 
even in this simple Keplerian picture, one may expect a 
dominance of the centrifugal force over gravity at suffi- 
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FIG. 20: Rotation rate /3b at the time of bounce for all models 
versus the precollapse central angular velocity f2 Cl i. The pro- 
genitor mass, the EoS, the precollapse differential rotation pa- 
rameter A, and the collapse dynamics are encoded as in Fig. 3. 
The lower horizontal line approximately separates pressure- 
dominated bounce models with and without strong prompt 
postbounce convection, while upper horizontal line marks the 
approximate transition between pressure-dominated bounce 
and centrifugal bounce. 

cicntly small r. In a more elaborate approach, employ- 
ing sequences of Newtonian self-gravitating equilibrium 
spheroids, Tohline [86] demonstrated that such a rota- 
tional barrier at which the collapsing core becomes cen- 
trifugally stabilized indeed exists in the context of stellar 
core collapse. This rotational barrier marks the hard up- 
per limit for the contraction of the inner core, hence also 
puts an upper limit /3 r b on the rotation rate that can 
be reached when varying Q c ,i for a given combination of 
precollapse degree of differential rotation and progenitor 
structure. 

Tohlinc's qualitative conclusions have been confirmed 
by multiple numerical studies of rotating collapse (see, 
e.g., [12-14, 16, 25] and our present work) while the quan- 
titative results, in particular the analytic critical rotation 
rate for centrifugal stabilization of collapse, do not hold 
for a dynamical collapse situation and must be deter- 
mined via nonlinear hydro dynamic simulations [18]. 

In Fig. 19 we plot the time evolution of the rotation 
rate (3 for a sequence of rotating collapse models with in- 
creasing precollapse central angular velocity f2 Ci ; while all 
other model parameters are kept fixed. All models reach 
their maximum rotation rate /3 max close to the time of 
core bounce, hence /3 max — /3b- After bounce, the inner 
core re-expands and settles into a new quasi-equilibrium 
configuration with /3 p b < /%■ Slowly to moderately 
rapidly rotating models experience little rotational sup- 
port, and in those cases /3b increases roughly linearly with 
51 c .i (see also Table III). For higher values of fi Cj i, cen- 
trifugal forces become relevant and /3b saturates at /3 r b as 
the models start to bounce centrifugally. For the s20A3 



FIG. 21: Rotation rate /3 p b in the late postbounce phase for 
all models versus the precollapse central angular velocity f2 c ,i. 
The progenitor mass, the EoS, the precollapse differential ro- 
tation parameter A, and the collapse dynamics are encoded 
as in Fig. 3. As in Fig. 20, the horizontal lines again ap- 
proximately mark the boundaries between different bounce 
dynamics. 



sequence with the Shen EoS considered here we deter- 
mine /3 r b to be ~ 23%. 

Figs. 20 and 21 provide an overview of the dependences 
of /3b and /3 p b, respectively, on f2 Cj ; for our entire model 
set as listed in Table III. Models that start out in es- 
sentially solid-body rotation (Al) never reach a /3b in 
excess of ~ 10% (with the maximum obtained in model 
s40AlO13). With increasing f2 Ci ; such rigidly-rotating 
cores become eventually fully centrifugally supported al- 
ready at the onset of collapse and do not collapse at all. 
Differentially rotating models may have higher values of 
f2 Ci i and thus a more rapidly rotating center, while the 
core is still allowed to collapse. As the collapse proceeds, 
electron capture reduces the pressure support and the 
size of the homologously collapsing inner core stays suffi- 
ciently small that centrifugal forces can become dynam- 
ically relevant only in the final phase of collapse (see the 
discussion in Section IV B). Thus, for our model set, 
the most differentially rotating configuration A3 leads to 
the highest values for /3b and /3 p b- A centrifugal bounce 
near the rotational barrier occurs only in a small sub- 
set of very rapidly (fi Cj i > 6.5 rad s _1 ) and differentially 
(A2/A3) rotating models, generally at /3b > 20.5%. 

At a fixed precollapse degree of differential rotation 
and f2 Cj i, /3b and /3 p b increase with a more massive and 
radially extended progenitor iron core (cf. Table I). This 
is analogous to the systematics found for the rotational 
enhancement of the inner core mass M- lc b at bounce (see 
Fig. 4). 

The dependence of both /3b and /3 p b on the EoS is 
small and shows little systematic trend. The Shen EoS, 
on the one hand, systematically yields a more massive 
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FIG. 22: Peak value |/i| max of the gravitational wave ampli- 
tude at 10 kpc distance for the burst signal (neglecting possi- 
bly larger contributions from postbounce convection at later 
times) for all models versus the rotation rate /% at the time 
of bounce. At slow to moderately rapid rotation, |/i| max is 
proportional to /3b to high accuracy (as marked by the dotted 
line with a slope of 1 in the log-log plot), while for /3b > 10% 
centrifugal effects reduce |/i| max . 



and more extended inner core that bounces with more 
dynamically relevant angular momentum than one ob- 
tained with the LS EoS. The LS EoS, on the other hand, 
leads to more compact configurations, which provide for 
stronger centrifugal spin- up in the final phase of collapse. 
The competition between these two effects results in the 
nonsystematic difference between the two EoSs seen in 
Figs. 20 and 21. 

The centrifugal barrier is also evident in Fig. 22, where 
we plot the dependence of the peak value |/i| max of 
the gravitational wave burst against the rotation rate 
/3b at bounce. It is noteworthy that centrifugal effects 
are responsible for an upper limit in |/i| max even be- 
fore the maximum rotation rate /3b ~ 25% is reached, 
which reflects the observation that the highest values of 
Hmax ~ 10 20 at 10 kpc distance are obtained for models 
which still undergo a pressure-dominated bounce, albeit 
at rapid rotation with /3b ~ 10%. Below these rotation 
rates, |/i| maat scales linearly with /3b with remarkable pre- 
cision, which is important information for the inversion 
problem in the case of a detection. We find a similar lin- 
ear dependence of |/i| ma x on the postbounce rotation rate 
/3 p b- In that case, however, the linear correlation is not as 
precise for low rotation rates (as /3 p b is rather sensitive 
to angular momentum redistribution due to convection 
after core bounce) and, in addition, the scaling becomes 
approximately quadratic well before |/i| max reaches its 
upper limit. 



B. The prospects for dynamical high-/3 instability 
in iron core collapse 

We find that none of our models surpass the threshold 
rotation rate /3d yn for the classical dynamical instability 
(see Table III). The overall largest /3 of ~ 25% is reached 
by model s40A3O15, which has the most massive and 
extended progenitor iron core (see Table I) in combina- 
tion with the strongest precollapse degree of differential 
rotation and highest precollapse central angular velocity 
considered in this study. This value of /3b ~ 25% comes 
close to the numerically obtained instability threshold of 
Pdyn > 25.5% reported in [77], but is maintained only for 
a very short time, since the core rebounds and settles at 
a more expanded quasi-equilibrium state after bounce. 
Accordingly, its postbounce rotation rate /3 p b is ~ 22%, 
and thus this model is unlikely to become subject to a 
dynamical high-/3 bar-mode instability. As portrayed by 
Fig. 21, the models with less extreme precollapse condi- 
tions in general reach a /3 p b significantly below ~ 20%. 

Based on the results from our extensive set of sim- 
ulations, we consider it unlikely that a proto-neutron 
star in nature develops a high-/3 dynamical instability 
at or early after core bounce. On the other hand, during 
its cooling to the final cold and condensed neutron star, 
the proto-neutron star contracts, and, if angular momen- 
tum is conserved and not redistributed or shed by other 
means (see, e.g., the discussion in [31, 59]), spins up on 
a timescale of seconds to minutes. While many of the 
proto-neutron stars in our model calculations could the- 
oretically reach /3d yn , it is, however, more likely that the 
secular instability driven by dissipation or gravitational 
radiation back-reaction, which in proto-neutron stars has 
a growth timescale on the order of 1 s [87], will set in 
first, completely diminishing the chances for dynamical 
high-/3 instability even in the most rapidly rotating proto- 
neutron stars. 

Finally, we point out that it is in principle possible to 
construct precollapse conditions that lead to /3b and /3 p b 
above /3d yn - This may be achieved by increasing signifi- 
cantly the precollapse degree of differential rotation and 
Q c _i above the values used in our most extreme models. 
However, such configurations (including already the rota- 
tional setup A3 in our models) are very unlikely to arise 
in evolution scenarios of single massive stars, since stellar 
evolution proceeds sufficiently slowly for redistribution of 
angular momentum towards solid-body rotation to occur 
on nuclear-burning timescales [31, 43, 88]. 



C. Differential rotation in the proto-neutron star 
and its relevance for the low-/3 dynamical instability 

The low-/3 dynamical instability appears to develop ex- 
clusively in differentially rotating fluid bodies and has 
been reported to occur even for rotation rates as low as 
~ 1%, provided the degree of differential rotation is suf- 
ficiently large [80]. 
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FIG. 23: Radial profile of the angular velocity f2 in the equa- 
torial plane at 20 ms after the time of core bounce for a rep- 
resentative subset of the models listed in Table III. Note that 
the inner core is in approximate solid body rotation out to 
about 10 km while the outer parts of the proto-neutron star 
and the postshock region rotate strongly differentially. The 
dotted lines mark the approximate range for the characteristic 
angular frequency fi c har- 

The nature of the low-/3 instability remains to be de- 
termined in detail, yet it has been suggested [81] that 
it is a type of dynamical shear instability that operates 
on the shear energy stored in differential rotation and 
radially redistributes angular momentum via the genera- 
tion of an azimuthal (nonaxisymmetric, spiral) structure 
that propagates outward in radius [25, 83]. In this pic- 
ture, nonaxisymmetric structure is generated by transfer 
of rotational energy from the axisymmetric background 
fluid to an azimuthal fluid mode at the location where 
the background angular velocity matches the mode pat- 
tern speed (i.e., at the corotation point). This proposed 
corotation mechanism suggests a close relationship of the 
low-/3 instability observed in simulations of stellar mod- 
els with dynamical instabilities in disks such as those 
described by Papaloizou and Pringle [89] . 

The importance of differential rotation for the low-/3 in- 
stability in stars can now be understood by the combina- 
tion of two important factors: First, differential rotation 
provides the reservoir of shear energy that can be tapped 
to generate the nonaxisymmetric structure. Second, de- 
spite a relatively low global rotation rate (3, differential 
rotation allows the central regions of a star to rotate suf- 
ficiently rapid to be in corotation with the lowest-order 
unstable modes that have pattern speeds of 0(27r/Td yn ), 
where 

t ^Hgm (24) 

is the dynamical timescale of the rotating star set by the 
Keplerian angular velocity [22, 79]. 

Since solid-body rotation is the state of lowest ro- 
tational energy, neutron stars are very likely to be- 
come rigidly rotating within at most a few dissipative 
timcscalcs during their post-supernova cooling evolution. 



Significant differential rotation may be expected in early 
merger remnants of binary neutron stars (e.g., [90]) and, 
importantly, is a consequence of rotating iron core col- 
lapse to a proto-neutron star investigated in the present 
work. 

In Fig. 23 we plot radial profiles of the angular velocity 

in the equatorial plane at 20 ms after core bounce for 
several of our models. As a result of quasi-homologous 
contraction, the near uniform precollapse rotational pro- 
file of the inner core is essentially frozen in during col- 
lapse [59]. In the outer core, however, the collapse pro- 
ceeds supersonically, resulting in differential rotation at 
equatorial radii > 10 km. In all models shown in Fig. 23, 
O declines by about two orders of magnitude in the ra- 
dial interval from 10 to 200 km, and roughly obeys a 
power-law with an exponent in the range of —1.2 to —1.4. 
Generally, a stronger degree of precollapse differential ro- 
tation leads to a steeper radial decline of after bounce. 
When increasing f2 Cj i while keeping the degree of prec- 
ollapse differential rotation fixed, the outer core regions 
experience more centrifugal support during collapse, re- 
sulting in a shallower postbounce slope for f2 (cf. model 
S20A2O15 in Fig. 23). 

In general, we find that the central angular velocity 
fl c after bounce increases monotonically with the prec- 
ollapse value f2 c ,i- For our models we obtain values for 
il c in the nascent proto-neutron star between about 2 
and 6 rad ms -1 , which corresponds to central rotation 
periods of about 1 to 3 ms. Assuming a mass range of 
the proto-neutron star of ~ 0.6 to 0.8 M Q for the models 
considered here (see Fig. 4) and a fiducial radius of the 
inner core at bounce of ~ 20 km, we obtain dynamical 
times of 1.7-2.0 ms, which yield characteristic angular 
frequencies of f2 c har ~ 3-4 rad ms -1 . Since the lowest 
order unstable mode is likely to have a pattern speed of 
the order of f2 c har, most models whose angular velocity 
we plot in Fig. 4 may indeed have corotation points with 
an unstable mode, hence could undergo a corotation-type 
low-/? instability. Slow rotators (with f2 Ci j < 2 rad s _1 ) 
do not appear to reach a sufficiently high angular velocity 
in the inner proto-neutron star core to have corotation 
points with potentially unstable modes in the first sev- 
eral tens of milliseconds after bounce. However, this may 
change at later times when the proto-neutron star con- 
tracts and spins up. 

Finally, we point out that our discussion is based on a 
very rough estimate of the pattern speed for the lowest- 
order unstable azimuthal mode. More reliable estimates 
can be made via multi-dimensional pcrturbative analysis 
(see, e.g., [83] in the context of idealized models) or by 
performing a large set of numerical simulations in three 
dimensions, which we plan to carry out in a future study. 



VIII. SUMMARY AND CONCLUSIONS 

In this article we have presented results from a com- 
prehensive set of collapse simulations of rotating stellar 
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iron cores to proto-neutron stars, using the axisymmetric 
general relativistic hydrodynamics code CoCoNuT. Our 
simulations treat all the relevant physics of the collapse 
phase to good approximation. They include precollapse 
iron core profiles from stellar evolutionary calculations, 
a highly efficient approximate treatment of deleptoniza- 
tion, a microphysical finite-temperature EoS, as well as 
neutrino pressure contributions. Magnetic fields are not 
included, since their relevance in the collapse and early 
postbounce phases is very likely negligible in cores with 
realistic precollapse fields [5, 31, 33, 37]. 

The focus of our study is on procuring accurate and 
reliable waveforms of the gravitational wave burst sig- 
nal associated with core bounce and on understanding 
the dependence of the signal characteristics on progeni- 
tor star mass, precollapse rotational setup, and nuclear 
EoS. To this end we have performed the to-date most 
extensive parameter study of this scenario, covering with 
more than 100 model calculations the parameter space 
spanned by (1) progenitor mass and model profile (zero- 
age main sequence masses from 11.2 to 40 Mq, presuper- 
nova models with and without rotation), (2) rotational 
configuration (slow and uniform to rapid and differential 
rotation), and (3) nuclear EoS prescription (from rela- 
tively soft to relatively stiff) . Importantly, the parameter 
space encompasses and even goes beyond all precollapse 
rotational configurations that are deemed realistic in the 
context of collapsing massive stars. 

A central result of this work is the finding that the 
gravitational wave burst from core bounce exhibits a 
generic waveform shape known as type I in the liter- 
ature [12, 13], independent of the model parameters. 
The multiple centrifugal bounce dynamics and the corre- 
sponding type II waveform found in previous, technically 
less complete studies (see, e.g., [12-14, 16]) do not occur 
in our models. 

We have demonstrated that all models with precol- 
lapse core angular velocities fl c ^ below ~ 5 rad s _1 (cor- 
responding periods longer than about 1 s) reach nuclear 
densities and experience a core bounce predominantly 
due to nuclear pressure effects. More rapidly rotating 
cores develop sufficient rotational support during collapse 
to undergo either a mixture of centrifugal and pressure- 
dominated bounce or a single centrifugal bounce at sub- 
nuclear densities. Centrifugal hang-up much below nu- 
clear density or multiple, damped harmonic oscillator- 
like centrifugal bounces do not occur. Therefore, these 
models also exhibit a type I waveform. The detailed 
analysis of the collapse dynamics presented in this pa- 
per reveals that the combined effects of general relativity 
and deleptonization lead to an increased destabilization 
of the collapsing core, result in a relatively small radius 
and mass M; c .b of the sonically-connected inner core at 
bounce (but not small enough to show the type III wave- 
form associated with rapid collapse found in some previ- 
ous simplistic models), and diminish the dynamical im- 
portance of centrifugal forces during collapse. 

The key parameter which determines the peak ampli- 



tude | h | max of the gravitational wave burst has turned 
out to be the precollapse central angular velocity f2 C; i. 
Slowly rotating cores with fi C) i < 1 rad s _1 produce fee- 
ble peak amplitudes on the order of 10 -22 at a distance 
of 10 kpc. More rapidly rotating cores with 1 rad s _1 < 
^c,i ^5 6 rad s _1 develop stronger quadrupole deforma- 
tions and have a rotationally-increased mass M; c .b at 
bounce, resulting in sizeable peak amplitudes in the range 
of 5 x 10~ 22 < |/i| max < 10~ 20 . The peaks of the wave- 
form spectrum of such cores cluster in frequency space in 
the interval of 650-800 Hz. At larger f2 c .i, centrifugal ef- 
fects become strong, significantly decelerate collapse and 
bounce, and even lead to a purely centrifugal bounce in 
a subset of models. This results in a general decrease 
of | h | m ax and a shift of the waveform's spectral peak to 
frequencies below ~ 400 Hz at high f2 Cj i. 

We have also shown that, in addition to f2 Ci j, the prec- 
ollapse core mass in combination with the electron frac- 
tion sets the mass Mj Cj b of the inner core at bounce, is an 
important quantity influencing the strength of the grav- 
itational burst. Since more massive progenitors gener- 
ally (though with notable non-monotonicity in the mass 
range from about 18 to 23 Mq) form larger iron cores, 
we observe in our model series a general trend to bigger 
Mic,h and larger |/i| max with increasing progenitor mass 
if all other parameters are kept constant. For instance, 
the 40 Mq progenitor yields values of |/i| max which are 
up to 4 times larger than for the lower-entropy 11.2 Mq 
counterpart with the same rotational configuration. 

The variations in the degree of differential rotation con- 
sidered in this study have only a minor impact on the col- 
lapse dynamics and burst waveform amplitude. Increas- 
ing differential rotation at fixed O c i generally lowers the 
centrifugal support of outer core regions. However, since 
the dynamically most relevant inner core at bounce con- 
sists of only ~ 0.5 -1M Q located within about 1000 km 
at the onset of collapse, the effects of differential rotation 
on the gravitational wave burst are small. 

Our results further indicate that the nuclear EoS has 
little influence on the gravitational wave burst signal. For 
a given precollapse configuration, a softer nuclear EoS 
yields higher densities at bounce and postbounce times 
with shorter variation timescales of the quadrupole mo- 
ment, but also leads to greater inner core compactness. 
In our simulations, the two effects generally cancel, lead- 
ing to no systematic trend in the peak waveform am- 
plitude | h | max with the EoS. The peak of the waveform 
spectrum, however, shifts to higher frequencies in the 
case of a softer EoS. For the models considered here, this 
frequency shift amounts to ~ 5.5% on average for mod- 
els undergoing pressure-dominated bounce. It is signifi- 
cantly smaller for models bouncing at subnuclear densi- 
ties under the influence of centrifugal effects. 

If situated within our Galaxy, a large fraction of our 
models are comfortably detectable by current gravita- 
tional wave detectors with a signal-to-noise ratio of up to 
6 in the most optimistic case (which is obtained for the 
most rapidly rotating models that still undergo pressure- 
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dominated core bounce). Advanced detectors could ob- 
serve them easily out to ~ 100 kpc and up to several 
10 Mpc for third-generation detectors. 

While such a gravitational wave signal may per se be 
detectable, the extraction of detailed physical informa- 
tion from the signal (i.e., solving the "inversion prob- 
lem" ) from the signal will be a formidable task. The very 
generic morphology of the burst waveforms and the clus- 
tering in frequency space of most models make it seem 
unlikely that a pure waveform-template-based inversion 
(as, e.g., carried out in [91] using the waveforms of [16]) 
can be successful for determining key physical parame- 
ters to significant precision. Our results, however, sug- 
gest that based on |/i| ma x and the peak frequency /max 
of the waveform spectrum alone, it should be possible to 
discriminate between purely pressure-dominated bounce 
(small to large |/i| ma x at frequencies / max significantly 
above 500 Hz) and centrifugal bounce (large |ft| max at 
frequencies / max significantly below 500 Hz). Further- 
more, we find that for not too rapid rotation |/i| max can 
be directly used to extract the rotation rate /?b at bounce 
to good precision. 

Making use of the extensive set of postbouncc ro- 
tational configurations obtained with our simulations, 
we have also studied the prospects for the development 
of nonaxisymmetric rotational instabilities in nascent 
proto-neutron stars. We find that the rotational bar- 
rier imposed by centrifugal forces prohibits the spin-up 
to rotation rates necessary for the classical dynamical 
bar-mode instability at high values of j3. We find, how- 
ever, that a large subset of our postbouncc models ex- 
hibits sufficiently differential and rapid rotation to be- 
come subject to the recently discovered low-/? instability. 
Still, three-dimensional simulations as in [17, 19, 22, 38] 
will be necessary to provide conclusive tests of our predic- 
tions. Furthermore, the interaction and competition of 
the low-/? instability and other instabilities operating on 
the shear energy of differential rotation, for instance the 
magneto-rotational instability (see, e.g., [33, 92]), remain 
to be studied. 

Finally, we point out that this study may be regarded 
as part - with the presently highest level of sophistication 
- of a multi-decade effort of our groups [11-14, 16-18] to 
provide reliable estimates for the gravitational wave burst 
emission associated with rotating core collapse and core 
bounce. The waveforms presented here are for the first 
time not only accurate (i.e., numerically converged), but 
reliable and robust, since our calculations take into ac- 



count all the necessary physics, including general relativ- 
ity, deleptonization and a microphysical EoS. All wave- 
forms are available for download in various formats in a 
publicly accessible waveform catalog [42]. 

We point out that the gravitational wave emission 
process considered in this work operates at measurable 
strength only if the progenitor core is rotating a lot more 
rapidly than expected for ordinary iron cores (see, e.g., 
[31, 59]). In slowly rotating core-collapse supernovac, 
turbulent convective overturn, instabilities of the accre- 
tion shock, and, possibly, proto-neutron star pulsations 
are likely to be the dominant emission processes of gravi- 
tational waves. The characteristics of these emission pro- 
cesses are not as well understood and will require more 
extensive and precise modeling to provide accurate es- 
timates of the complete gravitational wave signature of 
core-collapse supernovac. 
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